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Chapter 1 

Foreword 


Sextractor for dummies is not a book in the famous series, I just liked the title. Sextraetor is not a 
toy for grownups, it is a extremely usefull and versatile astronomieal (and perhaps other disiplines) 
softwaretool. 

In the eourse of my PhD I learned a great deal about it by trial and error and from Ed Smith and 
many others. To have somewhere to explain every knob and dial is the motivation behind this man¬ 
ual.^ I hope to make it as aeeurate and eomplete and possibly even eomplete as possible. I hope it 
helps with whatever projeet you are working on. 

There has been quite some time between updates of this user manual. I apologize but eertain 
external faetors insisted I write a PhD thesis first. In the meantime, both the offieial doeumentation 
and the program have been updated. I hope I’ve now eaught up suffieiently that this is again a 
useful manual for the beginning and advaneed souree extraetor user. ^ 


* Yes my motto is “Never read a manual as long there are still buttons and levers left...” but in sextractor’s case you 
might want to make an exception. Besides I’m perfectly happy if you just skim it. 

^And if you find a problem, inconsistency or spelling erro...please drop me a line at 
holwerda® stsci.edu 
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Chapter 2 

Introduction 


Figure 2,1: ??Just playing around with mkobject in iraf, source extractor and photoshop. 

Source Extractor ^ is used for the automated detection and photometry of sources in fits image- 
files. SE works on scans of phofographic plafes as well as CCDs. However I am assuming you are 
working wifh ccd dafa. 


If you goal is fo gef cafalogs of all fhe defecfed objecfs wifh reasonably good 
phofomefry from a FITS file of processed asfronomical imaging dafa, fhen SE 
is fhe insfrumenf you need.^ 


If you wanf REAEEY good phofomefry on only a few objecfs and you know where fhey are, 
do nof use SE. I’ll discuss how SE works and subsequenfly how fo sfarf if and which paramefers do 
whaf and whaf fhe differenf paramefers in fhe cafalogs mean. Some, buf by no means all, sfrafegies 
for ifs use are discussed and fo illusfrafe some examples from my own experiences will be given. 

This sfarfed as a bunch of notes of mine and Ed Smifh and is evolving slowly beyond fhaf. I 
hope fo provide for slighfly more insighf info SE fhan fhe official manual for fwo reasons: firsl, all 
fhe paramefers (bofh inpuf and oufpuf) will be discussed in one documenf, making if slighfly more 
complefe fhan eifher of fhe manuals and secondly if has DON’T PANIC wriffen in large friendly 
leffers on fhe cover. 


'Sometimes abbreviated to ’sex’, like the executable. If this funny to you a second time you might want to get out 


more. 






Chapter 3 

Pros and Cons of SE 


As I tried to point out in the introduetion and will emphasize oeeasionally SE is not suited for 
every astronomy projeet that needs photometry of objeets in a field. As everything, it has strong 
and weak points. 

The pros of SE are listed in the manual but the most important ones are: 

1. Speed. SE is made to go through data quiekly. And if you’re trying to beaver through several 
square degrees of data, speed is GOOD. 

2. The eapaeity to handle large fits files. SE is eoded up so fhaf if’ll lake if a pieee al a lime. 
Again good for Ihe sky-eafers among us. 

3. Works on CCD and seanned phofographie plale dafa. Niee if you happen lo have Ibis kind 
of dafa. 

4. Does deeenl pholomelry. 

5. Robusl, il’ll run wilh idiolie inpul. 

6 . Conlrollable, mosl steps ean be influeneed by user. 

7. The possibility lo aeeepf user speeified flag images or weighl images. 

8 . Oulpul paramelers and Ihe order in which Ihese are lisled, are specified by Ihe user. 

8 . Oulpul image-files depicting aperfures, deleclions and more. 

9. The possibility lo delecl sources in one image and do Ihe pholomelry in anolher. 

10. There is follow-up software lo decomposition of galaxy profiles. Source identification and 
a lirsl selection can be done wilh source exlraclor and eilher GIM2D or GAEEIT can Ihen 
follow-up. 
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However SE has some drawbacks. It was made for speedy use and in some case accuracy has 
been sacrificed for speed on purpose. So here is the other side of the coin: 

1. Only as good as its settings. SE is dependent on some of it’s setting and these are crucial for 
the detection and photometry. It will run on just about any set of input parameters but give 
back output that may be total bogus. 

2. Manuals are outdated and incomplete. This handbook is written as a remedy for that but by 
a user and not the person who wrote the code. ^ 

3. Eimited accuracy. You’ll see this with the geometrical output parameters. These are com¬ 
puted (from moments), NOT fitted (which would be more accurate). 

4. Classification of objects is of very limited use. 

5. Breaks down in crowded fields eventually. 

6 . Corrections of photometry for the ’wings’ of object profiles is very rudimentary. 

But as long as you, the user are aware of these little drawbacks and use the SE input as a fist 
start to fit with GIM2D or something then all will be well. 


* So I’m NOT claiming completeness or correctness. 



Chapter 4 


How to install SE 


4.1 Installing version 2.2 

First you get the most recent version of SE from http://terapix.iap.fr/soft/sextractor/index.html 
Then you unzip and tar the file (UNIX): 


gzip -dc sex_2.2.2.tar.gz | tar xv 

That should leave you with a directory in the directory where you did this called sextrac- 
tor2.2.2/ with instructions on how to install in ’INSTALL’. Basically you go to the ’sextrac- 
tor2.2.2/source’ directory and type: 


make SEXMACHINE= 'machine type' 


where the ’machine type’ can be any of the following possibilities: 


aix 

{for 

alpha 

{for 

hpux 

{for 

linuxpc 

{for 

linuxp2 

{for 

linuxkV 

{for 

sgi 

{for 

Solaris 

{for 

SunOS 

{for 

ultrix 

{for 


IBMs RS6000 running AIX) 
DEC-ALPHAs with Digital UNIX) 
HP/UX systems) 

PCs running LINUX, using gcc) 
Pentium2/3/4 PCs running LINUX) 
Athlon PCs running LINUX) 

SGI platforms) 

SUN-Solaris machines) 

SUN-OS platforms) 

DEC stations running ULTRIX) 


The SE manual is available in postscript format to you in the ’sextractor2.2.2/doc’ directory. 

Congrats, you have Source Extractor available as an executable ’sex’ in the ’sextractor2.2.2/source’ 
directory. The ’make’ file fries fo make a shorfcuf fo fhis execufable in your ’home’ directory. If fhis 
fails, fry making an alias for fhe command or simply fype fhe whole path/wherever/sextractor2.2.2/source/sex 
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4.2 Installing version 2.3 and up 

The simplest way to compile this package is: 


CHAPTER 4: How TO INSTALL SE 


1 . ‘cd’ to the directory containing the package’s source code and type ‘./configure’ to configure 
fhe package for your sysfem. 

2. Type ‘make’ fo compile fhe package. 

3. Type ‘make insfalT fo insfall fhe programs and any dafa files and documenfafion. 

Buf if you use linux, fhe rpm files are also available. 

4.3 Scisofl 

Af fhe European Soufhern Observafory, a package of all inferesfing scienlific soflware is being kepf 
up fo dale and easy fo insfall. Source exlraclor is pari of Ihis package and since fhe inslallafion is 
almosl plug-and-do-science (for fhe mac al leasl...). I highly recommend if, especially fhe Mac.^ 
gel if here: 

\protect\vrule widthOpt\protect\href{http://www.eso.org/science/scisoft/}{htt 
\protect\vrule widthOptXprotect\href{http://www.stecf.org/macosxscisoft/}{htt 


*Ok....rve plugged this enough. 



Chapter 5 

How SE works 


The source extractor package works in a series of steps. It determines the background and whether 
pixels belong to background or objects. Then it splits up the area that is not background into 
separate objects and determines the properties of each object, writing them to a catalog. 

The background determination is treated in the official manual and in section FOI All the pixels 
above a certain threshold are taken to belong to an object. If there is a saddle point in the inten¬ 
sity distribution (there are two peaks in the light distribution distinct enough), the object is split in 
different entries in the catalogs. Photometry is done on these by dividing up the intensity of the 
shared pixels. There is an option to ’’clean ” the catalog in order to eliminate artifacts caused by 
bright objects. Afterward, there is a list of objects with a series of parameters measured (ellipticity, 
size etc.). These are classified info sfars and galaxies (everyfhing non-sfar) by a neural nefwork. 

The firsl sfeps are confrolled by a number of paramefers. How fo esfimafe fhe fhreshold? How 
much confrasf should fhere be fo splif an objecf? However, fhe classificafion by fhe neural nef- 
work only depends indirecfly on fhe paramefers confrolling fhe firsf sfeps. As fhe nefwork has been 
frained on ground based dafa, fhere mighf be some doubfs on fhe reliabilify of fhis classificafion as 
one swifches fo ofher passbands or insfrumenfs. So check up on fhis classificafion in fhe case of 
fainf or blended objecfs. 

Some of fhe sfeps in SE have maps associafed wifh fhem and fhese can be wriffen fo fifs files, 
fhe ’check’ images. 

Aside from fhe inpuf image, SE can handle weighfimages and flag-images fo mark fhe relative 
imporfance of pixels or fo flag bad ones. 

The filling of fhe Poinl-Spread-Euncfion is nol yef operalional bul if was already indicaled in 
fhe second flow diagram in manual version 2.1.3. 

Sfeps of SE: 

1. Measure fhe background and ils RMS noise (background and RMS maps). (6.2) 

2. Subfracf background. 
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3. Filter (convolve with specified profile). (6.3.2) 

4. Find objecfs (fhresholding). (6.3.1) 

5. Deblend defecfions (break up defecfion info differenl objecfs. (6.3.3) 

6 . Measure shapes and positions. 

7. Clean (reconsider defections, accounfing for confribufions from neighbors) (6.3.4) 

8 . Perform phofomefry. 

9. Classify/index level of fuzziness more sfar-like or galactic? 

10. Oufpuf Cafalog and ’Check’ Images 

Figure 5.1: SE flow diagram from fhe firsl manual. This shows fhe principal steps done by SE. 

Figure 5.2: SE flow diagram from fhe second manual. This shows clearly fhe many exfra options 
over a simple run of SE; fhe possibilities of weighl maps, flags, crosscorellafion cafalogs and many 
more. 




Chapter 6 


Using SE 


Source Extractor can be used in basically three ways: on a single file for both detection and pho¬ 
tometry, on two files, one for defecfion, fhe ofher for phofomefry and on fwo files wifh cross- 
idenfificafion in fhe cafalogs. 


6 .1 Using SE on one image 

SE needs a series of paramefers in order fo run and fhese can be given af fhe command line or in a 
configurafion file. So fo run SE on a single file wifh all fhe necessary paramefers in fhe configura¬ 
tion tile one fypes: 

sex image -c contigurafion_tile.fxf 

If fhere is no configuration tile given, SE will fry fo read ’defaulf.sex’ from fhe local direcfory. 
However, fhe paramefer values can be fed fo SE on fhe command line as well: 

sex image -c contigurafion_tile.fxf -PARAMETERl valuel -PARAMETER2 value2 

The names of fhe paramefers and fheir meanings and preferred values are discussed below. 


NOTE: if you use bofh a configurafion tile and command line paramefer inpuf, 
fhe command line inpuf fakes precedence over fhe configurafion tile value. 


6.2 Using SE on separafe images for defection and phofomefry 

This may sound weird buf fhe option in SE fo find all fhe objecfs in one image and fhen apply fhe 
aperfures and posifions found on anofher image can be quife useful. To use imagel for fhe defecfion 
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of sources and image2 for the photometry: 
sex image 1,image! -c configuration_file.txt 

Suppose you want different colors or color information from the series of images you have in 
different filters. It is in that case very convenient to use the same apertures. So provided the images 
are well aligned, ' the photometry done is essentially the same objects using the same apertures.^ 
The nice part is that the numbering in the catalogs in this dual mode and the numbering in the 
single mode on the first image are the same. It should spare you a lot of rooting around in the 
catalogs if you want to compare fluxes in differenl bands for insfance. This particular fealure is 
furlher discussed in fhe sfrafegy secfion. 

6.3 Crosscorrellafing cafalogs. (ASSOC) 

There is a fhird possibilify fo gel informalion on objecls lhal occur in separafe images, for insfance 
in overlapping fields in a survey or simulafed objecfs in fhe SE deleclions. Basically you run SE 
on one field and fake fhe X and Y posilions (in pixels) from fhe cafalogs and feed Ihem lo SE in a 
second run wifh a search radius and a priority (fhe brighlesl associalion or fhe nearesf? elc.) And 
fhe malches are printed in a new calalog. See fhe ASSOC parameters in fhe calalog configuration 
secfion. 

6.4 Mulli-exlenlion fils 

There is a new fealure in Source exlracfor. If now supporls mulliple extension fils files. EITS 
is now fhe ’slandard’ filelype^ for aslronomical observations and several images can be pari of a 
single fils-file. Many new dala-producls now include weighl and quality maps of fhe astronomical 
image 


'The bright sources are on the same pixels, check by loading both in saotng or another display package and then 
’blink’ between them. If you’re drizzling or stacking several exposures, use the same reference image. 

^If you want awfully good photometry then it might be good te realize that a point spread function correction is 
dependent on the filter used. But keep in mind that SE does not do THAT good photometry to start with. 

^As with standards, everyone seems to have their own. 

"^The Advanced Camera for Surveys on board HST for instance. 




Chapter 7 


SE Input: the Configuration File 


”On two occasions, I have been asked [by members of Parliament], ’Pray, Mr. Babbage, if you 
put into the machine wrong figures, will the right answers come out?’ I am not able to rightly 
apprehend the kind of confusion of ideas that could provoke such a question.” 

- Charles Babbage ^ 

As stated in the previous section, SE tries to read the configuration file default.sex or a file 
can be given on fhe command line.^ The default.sex can be found in fhe /sextractor-2.4.4/config/ 
directory. If gives a good sef of defaulls for SE fo use. 

The configuration file is good way of remembering which paramefers you used in running SE 
and you do nof have to resef all fhe parameters when running SE over a series of files. The config- 
urafion file is an ASCII file (plain fexf) wifh fhe name of fhe paramefers and fhe value on separafe 
lines. A commenf line begins wifh ’#’ and ends wifh fhe end-of-line. 

The parameters are lisfed alphabefically in fhe manual. I discuss fhem here in fopical order. 
Inpul paramefers for SE can be roughly divided info fhe following categories: image information, 
background esfimafion, defection, pholomelry, cafalogs and SE running paramefers. 


7.1 Image Information 

SE gefs fhe posifional informalion from fhe EITS header buf mosf of fhe following paramefers musl 
be specified. GAIN is fhe rafio of fhe number of elecfrons to fhe number of ADU. The GAIN is 
dependenf on fhe fype of CCD you’re using and fhe insfrumenf in fronl of if (for inslance fhe WE2 
chip on fhe HST has a gain of 7). 

How fhe number you use here, fhe effeclive gain, relales to fhe insfrumenf gain is as follows: 


'Morale: input is inportant... 

^If the idea of having .sex files littering your hard disk is a little too randy for you, SE will read in any ASCII text 
files it is given as a configuration file. 
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Effeetive Gain 

Magnitude zeropoint 

Type of image 

gain X total exposure time 

zeropoint(l see) 

input image is e/s 

gain 

zeropoint(l see) - 1 - 2.5 
logio(exp. time) 

sum of N frames 

N X gain 

zeropoint(l see) - 1 - 2.5 
logio(av. exp. time) 

average of N frames 

2xNxgain/3 

zeropoint(l see) - 1 - 2.5 
logio(av. exp. time) 

median of N frames 


NOTE: There are different strategies you ean follow using the GAIN and the 
zeropoints, see the SE strategies seetion. 


The MAG_GAMMA is a relie of when this program was applied to photographie plates (seans 
thereof). SE v2.2.2 did not RUN, in my ease, without it speeified, even though it does not use it 
while proeessing your niee CCD data. Fortunately this seems fixed in a reeent update. If you have 
the misfortune that your data are still seanned in photographie plates, then this is the slope of the 
response funetion of the emulsion used on the plates in question. 


NOTE: SE (before v2.4.4) will not run without MAG_GAMMA speeified. 


What does MAG_GAMMA do? (from the offieial manual) 

Photographie photometry In DETECT TYPE PHOTO mode, SExtraetor as¬ 
sumes that the response of the deteetor, over the dynamie range of the image, 
is logarithmie. This is generally a good approximation for photographie den¬ 
sity on deep exposures. Photometrie proeedures deseribed above remain un- 
ehanged, exeept that for eaeh pixel we apply first the transformation 

I = /olof (7.1) 

where (7 = MAG_GAMMA is the eontrast index of the emulsion, D the original 
pixel value from the baekground-subtraeted image, and Iq is eomputed from the 
magnitude zero-point mo (speeified in ): 

One advantage of using a density-to-intensity transformation relative to the 
loeal sky baekground is that it eorreets (to some extent) large-seale inhomo¬ 
geneities in sensitivity. 
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The DETECT_TYPE specifies what type of data SE is handling, scanned photo plates or 
CCD data. Even with DETECT_TYPE set to CCD, older versions of SE will still need that 
MAG.GAMMA. 

MAG2ZEROPOINT is the zeropoint for the photometric measurements. This is again different 
if you use counts-per-second images as opposed to total counts images. The counts-per-second 
images have the zeropoint specified by the instrument handbook (depends on filter, instrument 
and type of ccd used). But in the case of a total counts image is the handbook value plus the 
2.5 logio(exposuretime). 

PIXEE_SCAEE is again something you hopefully know before you started to run SE. Eunny 
enough this is not read from the Eits header. So specify this! SE needs it only for the CEASS_STAR 
parameter (but still needs it). 


NOTE: New feature in SE (v2.4.4.). When this parameter is set to 0, SE uses 
the World coordinate fits information to compute the pixelscale. 


SATUR_EEVEE is the limit for SE to start extrapolating to get the photometry. However as 
soon as you hit something as saturated as that you might want to see if there is another way to 
determine the flux from that object. 

The SEEING_EWHM (Eull Width at Half Maximum) is important for that separation between 
stars and galaxies. Eike the PIXEE_SCAEE it is only used for the CEASS_STAR parameter.^ 


^Seeing is the blurring of the image as a result of atmospheric disturbances (cirrus clouds, turbulence etc.). An 
estimate of the seeing should be documented in either the header or the observation logs. If not, make something up, 
possibly inspired by that looks like a bright star. 
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Parameter 

Default 

Type 

Deseription 

FLAGJMAGE 

flag.fits 

strings 
(n < 4 ) 

File name(s) of the flag-image(s) . 

FFAG.TYPE 

OR 

keyword 

Combination expression for the 
method for flags on the same ob¬ 
ject: 



OR 

arithmetical OR, 



AND 

arithmetical AND, 



MIN 

minimum of all flag values. 



MAX 

maximum of all flag values. 



MOST 

most common flag value. 

GAIN 


float 

Gain (conversion factor in 
e”=ADU) used for error estimates 
of CCD magnitudes. 

DETECT.TYPE 

CCD 

keyword 

Type of device that produced the 
image: 



CCD 

linear detector like CCDs or NIC- 
MOS, 



PHOTO 

photographic scan. 

MAG.GAMMA 


float 

7 of the emulsion (slope of the re¬ 
sponse function). Takes effect in 

PHOTO mode only but NEEDS 
to be specified, even for CCD im¬ 
ages. 

MAG_ZEROPOINT 

” 

float 

Zero-point offset to be applied to 
magnitudes. 

PIXEF_SCAFE 


float 

Pixel size in arcsec. (for surface 
brightness parameters, FWHM 
and star/ galaxy separation only). 

SATUR_FEVEF 

” 

float 

Pixel value above which it is con¬ 
sidered saturated. 

SEEING_FWHM 


float 

FWHM of stellar images in arc- 
sec. This quantity is used only 
for the neural network star/galaxy 
separation as expressed in the 

CFASS_STAR output. 
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7.2 Background Estimation 

SE estimates the background of the image as well as the RMS noise in that background, mapping 
both SE subtracts the estimated background from the photometry and uses the RMS to estimate 
errors. So the background is important in the rest of the SE run. 

BACK_SIZE regulates the estimate. In an area of the BACK_SIZE, the mean and the a of the 
distribution of pixel values is computed. Then the most deviant values are discarded and median 
and standard deviation a are computed again. This is repeated until all the remaining pixel values 
are within mean ±3 cj. If a dropped with less than 20% per iteration, the field is considered not 
crowded. 

The value for the background in the area is: 

• the mean in the non-crowded case 

• 2.5 X median - 1.5 x mean in the crowded case 

Both the mean and the median are the ones computed in the last iteration. The mean is the 
average and the median is the average of all the values except the most extreme one. 

The background map is a bi-cubic-spline interpolation over all the area’s of size BACK_SIZE, after 
filtering. So obviously, the choice of BACK_SIZE is very important, too small and the background 
estimate will be partly object flux, too large and small scale variations cannot be taken into account. 
The effect of different BACK_SIZE values is illustrated in figure ITTI 


NOTE: The BACK_SIZE parameter determines the background map. Estimate 
the average size of the objects in pixels and make sure the BACK_SIZE is larger 
than that. 


But before the fit to the background values is done, the values can be smoothed: BACK_EIETERSIZE 
is the median filter for the background map. Effectively, you smooth the background image over 
this number of meshes to get rid of the deviations resulting from bright or extended objects. The 
effect of different BACK_EIETERSIZE values is illustrated in figure ??. 

BACK_TYPE is only used if you do not want SE to go off and estimate the background but use 
only one constant value supplied by you in BACK_VAEUE. 

The background computed above is first used in the detection of objects. Subsequently this 
background value can also be used for the photometry. But to get accurate background values 
for the photometry, the background can be recomputed in an area centered around the object in 
question. To recompute, set BACKPHOTO_TYPE to EOCAE and pick an BACKPHOTO_THICK 

''CHECKIMAGE.TYPE BACK GROUND and BACKGROUND_RMS if you want to inspect them. See also the 
checkimage section fsection l7.7.3> . 
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to match your tastes (generally speaking, somewhat larger than the objects in questions would be a 
good idea). 


NOTE: the RMS as determined from the BACKGROUND_RMS map will be 
used in more than just the photometry, the thresholds for detection and analysis 
can be dependent on it. 


NOTE: if you want to subtract the background and not have SE do this for you, 
set BACK_TYPE to MANUAE and BACK_VAEUE to 0.0,0.0 


Figure 7.1: The HDF-N with the background estimates using different BACK_SIZE values. 
Clearly a BACK_SIZE smaller than the biggest objects will result in a too much variation in the 
background estimate. A overlarge BACK_SIZE will not account for variations in background. 
Grayscales for the background field are the same. 
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Parameter 

Default 

Type 

Deseription 

BACK_SIZE 

- 

integers 

Size, or Width, Height (in 



(n < 2) 

pixels) of a baekground 
mesh. 

BACK_FIETERSIZE 

- 

integers 

Size, or Width, Height (in 



(n < 2) 

baekground meshes) of 
the baekground-filtering 
mask. 

BACK_TYPE 

AUTO 

keywords 

What baekground is sub- 



(n < 2) 

traeted from the images: 



AUTO 

The internal interpolated 
baekground-map. In the 
manual it says “INTER- 
NAE” here but the key¬ 
word is AUTO. 



MANUAE 

A user-supplied eonstant 
value provided in BACK 
VAEUE. 

BACK_VAEUE 

0 .0,0.0 

floats 

in BACK TYPE MAN¬ 



(n < 2) 

UAE mode, the eonstant 
value to be subtraeted 
from the images. 

BACKPHOTO.THICK 

24 

integer 

Thiekness (in pixels) of 
the baekground EOCAE 
annulus. 

BACKPHOTO.TYPE 

GEOBAE 

keyword 

Baekground used to eom- 
pute magnitudes: 



GEOBAE 

taken direetly from the 
baekground map. 



EOCAE 

reeomputed in a reetangu- 
lar annulus around the ob¬ 
ject. 


Figure 7.2: The HDF-N with the baekground estimates using different BACK_FILTERSIZE 
values. The mesh of baekground estimates (BACK_SIZE = 64) is smoothed over the 
BACK_FIETERSIZE before interpolating into the BACKGROUND map. 
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7.2.1 Weight Images 

The individual pixels in the deteetion image ean be given relative importanee by using a weight for 
eaeh of them. Different options are available: the baekground as determined by SE or an external 
weight image. If an external weight image is given, it has to be speeified what kind it is; a varianee 
map or a rms map or a weight map, from whieh a varianee map should be derived. 

The weight for eaeh pixel is derived as follows: 


weight =-=-; 

variance rms 


WEIGHT.TYPE MAP.WEIGHT is direetly taken as the weight, WEIGHT.TYPE MAP.VAR 
is inverted and WEIGHT_TYPE MAP_RMS is squared and inverted. 

Reasons for using a weight image are various but to give you an idea: SE ean ignore pieees of 
the image this way, use a flat-field in it’s photometry or edited results from a previous run. ^ 

Users of the stsdas drizzle paekage should eheek out the appendix on drizzle’s weight image. 

These are the parameters eontrolling it: 


I ean’t really fathom why WEIGTH_GAIN isn’t simply an option in WEIGHT_TYPE but it 
isn’t. So there. 


WEIGHT JMAGE is the input parameter where the fits file is given whieh is to be used as 
weight map of the type defined by the WEIGHT_TYPE parameter. 


With WEIGHT.TYPE set on BACKGROUND, the ’eheekimage’ (output image with CHECKIM- 
AGE_TYPE set on BACKGROUND) of a previous run ean be used. MAP_RMS ean be for instanee 
derived from known noise eharaeteristies of the instrument or given by other programs used in data- 
reduetion (for instanee the ’drizzle’ paekage for HST data). 


^Please note that I do not specify by what kind of black magic you obtained these weight images. They may be the 
result from a previous SE run but a ccd dark-frame might be used. Or perhaps your data-reduction scheme will produce 
a good estimate of the rms. To be used with caution! 
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Parameter Default Type Deseription 

WEIGHT_GAIN Y boolean If true, weight maps are eonsid- 

ered as gain maps. 

WEIGHT JMAGE weight.fits strings Eile name of the deteetion and 

(n < 2) measurement weightimage , re- 
speetively. 

WEIGHT_TYPE NONE keywords Weighting seheme (for single 

(n < 2) image, or deteetion and mea¬ 

surement images): 

NONE no weighting, 

BACKGROUND 

varianee-map derived from the 
image itself, 

MAP_RMS varianee-map derived from an 
external RMS-map, 

MAP_VAR external varianee-map, 
MAP.WEIGHT 

varianee-map derived from an 
external weight-map. 
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This is the description from the SE manual v2.3: 

NONE: No weighting is applied. The related WEIGHT IMAGE and 
WEIGHT THRESH (see below)parameters are ignored, (quick toggle so 
you can see if including the weight image has any effect) 

BACKGROUND: the science image itself is used to compute internally 
a variance map (the related WEIGHT IMAGE parameter is ignored). 
Robust (3 -clipped) variance estimates are first computed within the 
same background meshes as those described in x??12. The resulting 
low-resolution variance map is then bicubic-spline-interpolated on the 
y to produce the actual full-size variance map. A check-image with 
CHECKIMAGE TYPE MINIBACK RMS can be requested to examine 
the low-resolution variance map. 

MAP_RMS: the FITS image speci ed by the WEIGHT JMAGE le name 
must contain a weightmap in units of absolute standard deviations (in 
ADUs per pixel). 

MAP_VAR: the FITS image speci ed by the WEIGHT JMAGE le name 
must contain a weightmap in units of relative variance. A robust scaling 
to the appropriate absolute level is then performed by comparing this 
variance map to an internal, low-resolution, absolute variance map built 
from the science image itself. 

MAP.WEIGHT: the FITS image specified by fhe WEIGHT JMAGE file 
name musf confain a weighf-map in unifs of relafive weighfs. The dafa 
are converfed fo variance unifs (by definifion variance / l=weighf), and 
scaled as for MAP_VAR. MAP_WEIGHT is fhe mosf commonly used 
fype of weighf-map: a fiaf-field, for example, is generally a good ap- 
proximafion fo a perfecl weighf-map. 

In facl, if you wanf fo work wifh Weighl maps (since you fhink you have some- 
fhing fhaf mighf be used as a weigfh map), please read section 7 in fhe official 
manual. 
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7.3 Finding and Separating Objects 

SE considers every pixel above a certain threshold (to be specified by YOU ^ , directly or indirectly) 
to be part of an object. The ’deblending’ is the part where it figures ouf which pixels or parfs of 
pixels belong fo which objecfs. 

7.3.1 Defecfion; Thresholds 

The fhreshold paramefers indicate fhe level from which SE should sfarf freafing pixels as if fhey 
were pari of objecfs, delermining parameters from Ihem. There are Ihree requiremenls for a candi¬ 
date objecfs: 

1. All fhe pixels are above fhe DETECT_THRESH. 

2. All fhese pixels are adjacenf fo each ofher (eifher corners or sides in common). 

3. There are more fhan fhe minimum a number of pixels (specified in DETECTJVIINAREA). 

ANAEYSIS-THRESH is jusl fhe fhreshold for CEASS STAR and EWHM, all fhe ofher param- 
elers are determined from fhe DETECT_THRESH. 


NOTE: All fhe OTHER analysis (phofomefry and fhe like) is done wifh fhe 
DETECT.THRESH! 


They can bofh be specified in Ihree ways: 

• In Surface Brighfness: SBlimil,SBzeropoinl e.g. DETECT_THRESH 23.5, 24 

• In ADU’s (sef THRESH.TYPE ABSOEUTE) e.g. DETECT.THRESH 1.2 

• Relative fo background RMS (sef THRESH.TYPE REEATIVE) e.g. DETECT.THRESH 
1.2 (fhis is fhe initial sefup in fhe defaull.sex file.) 

The fhreshold in surface brighfness mu {mag/arcsec^ ) needs a calibration Zero-poinf {mag/arcsec^ 
corresponding fo 0 counls). Nole lhal fhis can/will be differenl from fhe MAG_ZEROPOINT value. 


NOTE: This whole SB fhreshold sluff seems quite fhe fhing unfil you realize 
SE jusl does fhis: 

_ ^^limit ~S^zeropoint 

thresh = 10 2.5 


®yes YOU! You did it to yourself, you and only you and that’s what really hurts. (Radiohead) 

^For the HST instruments, the zeropoints can be obtained using CALCPHOT in the stsdas package in IRAK Re¬ 
member that these are dependent on instrument and filter. 
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With THRESH_TYPE set to ABSOEUTE, the threshold is set to the same number of ADUs 
aeross the image. If THRESH_TYPE set to REEATIVE (the default), the the threshold is that 
number of baekground RMS standard deviations above the baekground value. This is niee and 
flexible but sensitive to the baekground estimation! 

The DETECT JMINAREA is the minimum number of pixels above the threshold required to be 
eonsidered an objeet. 


Figure 7.3: An objeet is defined as a series of bordering pixels above the threshold. Bordering 
pixels to the gray one are 1 through 8. Another pixels needs to be above the threshold and share a 
eorner or a side with another pixel in order to eonsidered part of the same objeet (unless they are 
deblended into different objeets). 
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Parameter 

Default 

Type 

Deseription 

THRESH.TYPE 

REEATIVE 

keywords 
(n < 2) 

REEATIVE 

ABSOEUTE 

Meaning of the DE¬ 
TECT THRESH and 
ANAEYSIS .THRESH 
parameters : 

sealing faetor to the 
baekground RMS, 
absolute level (in ADUs 
or in surfaee brightness). 

ANAEYSIS.THRESH 


floats (n < 
2) 

Threshold (in surfaee 
brightness) at whieh 
CEASS STAR and 
EWHM operate. 1 

argument: relative to 

Baekground RMS. 

2 arguments: mu 

(vnagj arcsed'’^ ), 

Zero-point (mag). 

DETECT.THRESH 


floats (n < 
2) 

Deteetion threshold. 1 
argument: (ADUs or 

relative to Baekground 
RMS, see THRESH 
TYPE). 2 arguments: 
R (mag.aresee” 2 ), 
Zero-point (mag). 

DETECTJVllNAREA 


integer 

Minimum number of 
pixels above threshold 
triggering deteetion. 


7.3.2 Filtering 

Before the deteetion of pixels above the threshold, there is the option of applying a filter. This filter 
essentially smooths the image. ^ 

There are some advantages to applying a filter before deteetion. It may help deteet faint, ex¬ 
tended objeets 

However it may not be so helpful if your data are very erowded. There are four types of filter to 
be found in the ./sextractor2.2.2/config/ directory, Gaussian, Mexiean hat, tophat and blokfunetion 
of various sizes, all normalized. 

This is what the helpful README in this direetory tells us: 

*The photometry is still being done on the original image, don’t worry. 
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Name 

default.conv 
gauss*.conv 

tophat*.conv 

mexhat*.conv 


block_3x3.conv 


Description 

a small pyramidal function (fast) 

a set of Gaussian functions, for seeing FWHMs between 1.5 
and 5 pixels (best for faint object detection), 
a set of ”top-hat” functions. Use them to detect extended, low- 
surface brightness objects, with a very low THRESHOLD, 
’’wavelets”, producing a passband-filtering of the image, tuned 
to seeing FWHMs between 1.5 and 5 pixels. Useful in very 
crowded star fields, or in the vicinity of a nebula. WARNING: 
may need a high THRESHOLD!! 

a small ’’block” function (for rebinned images like those of the 
DeNIS survey). 


The naming convention seems to be: name_seeingFWHM_size.conv. Both the seeing FWHM 
and the size are in pixels. So depending on what you are after, choose a filter and approximately 
your seeing. 


NOTE: filter choice and threshold choice are interdependent! 


Parameter 

Defaulf Type 

Description 

FILTER 

boolean 

If Irue, filtering is applied lo fhe 
dala before exlraclion. 

FILTER_NAME 

siring 

Name and palh of fhe file con- 
faining fhe filler definifion. 

FILTER_THRESH 

fioals (n < 
2) 

Lower and higher Ihresholds 
(in back-ground slandard devia- 
lions) for a pix-el lo be consid¬ 
ered in filtering (used for relina- 
fillering only). 


7.3.3 Deblending: separating into different objects 

Deblending is the part of SE where a decision is made whether or not a group of adjacent pixels 
above DETECT_THRESH is a single object or not. Suppose there is a little island of adjacent 
pixels above the threshold. It is an object or maybe several really close next to each other. So how 
does SE cut this up into different objects? First it defines a number of levels befween fhe fhreshold 
and fhe maximum counf in fhe objecf. This is sef by fhe DEBLEND_NTHRESH parameter. The 
levels are spaced exponenfially. 

SE fhen consfrucfs a ’free’ of fhe objecfs, branching every fime fhere are pixels above a fhresh- 
old separafed by pixels below if (see figure). A branch is considered a differenf objecf provided: 
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1. The number of eounts in the braneh (A in the figure) is above a eertain fraetion of the total 
eount in the entire ’island’. 

2. There is at least one other braneh (yep B!) above the same level that is also above this 
fraction. 

The fraction is defined in DEBLENDJVIINCONT. 

Figure 7.4: Illusfrafion of DEBLEND.NTHRESH and DEBLENDJVIINCONT. 

Paramefer Defaulf Type Description 

DEBLENDJVIINCONT - fioaf Minimum confrasf paramefer 

for deblending. 

DEBLEND_NTHRESH - infeger Number of deblending sub- 

fhresholds. 


7.3.4 Cleaning 

There is fhe option fo ’clean’ fhe lisf of objecfs of arlifacfs due fo brighf objecfs (sef CLEAN fo 
YES). All fhe defections are checked fo see if fhey would have been defecfed (i.e. exceeded fhe 
fhreshold efc.) if fheir neighbors were nol fhere. To do Ibis, fhe confribufions of fhe neighboring 
objecfs has fo be compufed. An esfimafe is made from a moffaf lighf profile. The Moffal profile is 
scaled and sfrefched fo til fhe neighbour’s profile. The confribufion fo fhe objecf from fhe wings of 
fhe Moffal profile is fhen subfracfed. 

The Moffal profile looks like Ibis: 


I{r) _ 1 

7(0) “ (1 + A: X r2)/3 

The CLEAN_PARAM is fhe /3 parameter in fhe above formula. 


NOTE: Decreasing CLEAN_PARAM yields brighter wings and more aggres¬ 
sive cleaning. 


The value for fhe CLEAN_PARAM should be belween 0.1 and 10. 


NOTE: In earlier versions of SE^*^, fhe Moffel profile was a Gaussian and fhe 
CLEAN_PARAM fhe slrelch factor for fhe EWHM. 

®Note that 0.01 is 1% of the flux when deflning this. 
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Cleaning would be more aggressive with a higher CLEAN_PARAM. This version of cleaning is 
explained in the original manual. 


Parameter 

Default Type 

Description 

CEEAN 

boolean 

If true, a cleaning of the catalog 



is done before being written to 



disk. 

CEEAN_PARAM 

float 

Efficiency of cleaning. 


7.4 Influencing Photometry 

After deblending the objects, SE performs astrometry (where stuff is), photometry (how bright stuff 
is) and geometric parameters (how stuff looks like). The astrometry cannot be influenced via input 
parameters (you just specify what kind of positions you want, easy). The photometry had a few 
input parameters associated with it; what to do with overlapping pixels, what is the zeropoint and 
how to apply apertures. To understand these, the AUTO photometry of SE has to be explained. The 
geometric parameters are mainly associated with the Kron (AUTO) photometry. 

The GAIN and MAG_ZEROPOINT, have been discussed at the image characteristics sec¬ 
tion. Of course they are needed for the photometry. The GAIN to convert counts to flux and the 
MAG_ZEROPOINT for the calibration of the magnitude scale. Also the BACKPHOTO_THICK 
and BACKPHOTO_TYPE give you influence on the way the background subtracted from the pho¬ 
tometry is estimated (see section IT^ . 

There are five different approaches in SB’s photometry; isophotal, isophotal-corrected, auto¬ 
matic, best estimate and aperture. ^' 


7.4.1 ISO 

In the above you defined above what threshold SE should do it’s photometry, with the estimated 
background as zeropoint. The pixels above this threshold constitute an isophotal area. The flux 
or magnitude determined from this (counts in pixels above threshold minus the background) is 
the isophot flux/magnitude. Apart from the threshold (DETECT_THRESH) and the background 
estimation, there is nothing to influence here. 

7.4.2 ISOCOR 

In real life however, objects rarely have all their flux within neat boundaries, some of the flux is in 
the “wings” of the profile. SE can do a crude correction for that, assuming a symmetric Gaussian 
profile for the object. This would be the isophot-corrected flux/magnitude. There is no parameter 
for you to influence this estimate. 

**SE v2.4.4 introduces a sixth, an ’automatic’ photmetry with the Petrosian radius. 
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How ISOCOR works (from the v2.3 manual) 

Corrected isophotal magnitudes (MAG ISOCOR) is a quick-and-dirty way for 
retrieving the fraction of flux in the wings of the object missed in the isophotal 
magnitudes. The latest version of the SE manual (v2.3) does not rate this as 
a very good correction, these values have been kept as an output option for 
compatibility with SE v2.x and SE vl^^. The assumtion is that the profiles 
of the objects are Gaussian. The fraction of the total flux enclosed within a 
particular isophote reads compared to the total flux^^: 

1 , / ^ 

1 _ _ in{l -r,) = — (13) 

I? 2 iso 

where A is the area and t the threshold related to this isophote. The 
relation can not be inverted analytically, but a good approximation to 
(error < 10^ for r] > 0.4) can be done with the second-order polynomial of 
t: 

At / A t\^ 

1-0.1961^-0.7512 — (7.4) 

-^iso \ ^iso ) 

A total magnitude mtot estimate is then mtot = niiso + 2.51og(7y) Clearly this 
first correction works best with stars (and maybe starclusters?) as these are 
the most Gaussian-like, ok with disk galaxies and not-so-ok with ellipticals as 
these have broader wings. 


Eixed-aperture magnitudes (MAG_APER) estimate the flux above the background within a circular 
aperture. The diameter of the aperture in pixels 

(PHOTOM_APERTURES) is supplied by the user (in fact it does not need to be an integer since 
each normal pixel is subdivided in 5 x 5 sub-pixels before measuring the flux within the aperture). 
If MAG_APER is provided as a vector MAG_APER[n], at least n apertures must be specified with 
PHOTOM_APERTURES. 


7.4.3 AUTO 


SE uses a flexible elli ptical apertu re around every detected object and measures all the flux in¬ 
side that, described in Kron ( 1980h . There are two parameters regulating the elliptical apertures: 


PHOT_AUTOPARAMS and PHOT_AUTOAPERS. 
The characteristic radius for the ellipse is: 


Sr/(r) 

S/(r) 


(7.5) 


Also known as the Kron radius (see section 18.5.11) . Erom the objects second order moments, the 
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ellipticity rj and position angle 9 are computed. The major and minor axes of the elliptical aperture 
are computed to be kri/e and ekri respectively. PHOT_AUTOPARAMS influences directly the 
estimation. The first is k factor mentioned above and the second is the minimum radius for a Kron 
ellipse. 

PHOT_AUTOAPERS are the minimum aperture diameters for the Kron photometry, estimation and 
measurement. These are circular! These are used in case the radius of the Kron aperture goes below 
the Rmin specified in PHOT_AUTOPARAMS. The values in the default.sex are probably best for 
most applications. A good way to check the setting on PHOT_AUTOPARAMS is to generate 
an ’APERTURES’ check-image and see if the Kron apertures overlap too much with those of 
neighboring objects. 

Parameter Default Type Description 

PHOT_AUTOPARAMS ^ floats (n = 2) MAG AUTO controls: 

scaling parameter k of 
the 1st order moment, 
and minimum Rmin (in 
units of A and B). 

PHOT_AUTOAPERS 0.0,0.0 floats (n = 2) MAG AUTO min¬ 
imum (circular) 

aperture diameters: 
estimation disk, and 
measurement disk. 


Neat trick: The Kron radius was introduced as a accurate flexible aperture 
that would capture most of the flux from an object. There is relation be- 
tween the Kron flux-meas urement and the total flux from an object given in 


Graham and Driven (120050 depending on the Sersic profile of the galaxy. See 


for the enclosed fraction table 18.21 


NOTE: In the 2.4.4 version of source extractor provides an alternative very sim¬ 
ilar to the AUTO photometry, the PETRO photometry parameters, determined 
within an aperture defined by the Petrosian radius. 


7.4.4 BEST 

With all this flexibility, you’d expect the Kron or AUTO photometry to be the best. However it can 
be influenced by nearby sources. Therefore there is a fourth option, MAG_BEST. This is usually 
equal to AUTO photometry but if the contribution of other sources exceeds 10%, it is ISOCOR. 

*"^Note that the elliptical apertures in figure ?? (Source Extractor written in sources? remember? Did you skim over 
it? you bastard!) are actual Kron apertures for the ’objects’. I cheated however and set PHOT.AUTOPARAMS to 1.0 
so the Kron radii would not be so big. 
















Typical Radii 


37 


Figure 7.5: An example of the different choices of the Kron factor (k) in PHOT_AUTOPARAMS. 
As the Kron aperture grows, it encloses more flux from other objects. Note how the checkimage 
set to APERTURES is a good check against this. Dodgy Kron apertures are dashed. 


NOTE: the characterization ’BEST’ is a bit misleading. I was advised from 
several sides NOT to use this photometry. Both AUTO and ISOCOR are at 
least consistent. 


7.4.5 Circular Apertures 

Then there is the user specified apertures. The aperture se diameters are specified in pixels in the 
PHOT.APERTURES. The fluxes, magnitudes and errors in these are specified in EEUXu^PER,MAGu^PER, 
EEUXERRu^PER and MAGERRu^PER. Please note that upon inclusion of N apertures, these pa¬ 
rameters are N columns each in the outputcatalog, corresponding to each specified aperture. 

The calibration documentation for the instrument used should give you a photometric correc¬ 
tion for a certain circular aperture and a point source. It tells you the relation between the flux 
captured within the aperture and all the flux from a star. Since part of the flux has been smeared out 
of the aperture by the PSP, you need to correct for this. If you have an extended object, the game is 
different. Portunately, there are other apertures^^ for this and even a correction for those to account 
for all the flux (but that also depends on the light profile of the galaxy). 

One of the geometric output-parameters is the half-light radius. The fraction of total light 
within this radius is specified in the PHOTJEUXPRAC parameter. 


Figure 7.6: Illustration of the different apertures possible: ISO, ISOCOR, AUTO and APER (user 
specified in PHOTu^^PERTURES 

7.5 Typical Radii 

There are several output option to describe the typical size of an object. Often they also define an 
aperture that is meant to capture an extended object like a galaxy. These are the input parameters 
that are important to them. The PEUXJIADIUS has been an output parameter for some time, the 
KRON and PETROSIAN radii appeared in version 2.4.4 of SE. 

7.5.1 Effective Radii 

SE has the option to put out radii containing a certain fraction of the light, (outputparameter 
EEUXJIADIUS) The default is 0.5 (the half light radius). PHOTJEUXERAC 0.2,0.5,0.9 will 
give three radii containing 20%, 50% and 90% of the light respectively. The effective radius output 
is discussed in section IS! 5. 31 


*^The Kron and Petrosian ones...see sections rS.S.II and ??. 
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Parameter Default Type Description 

PHOT_FLUXFRAC 0^ floats (n < Fraction of FFUX 

32) AUTO defining each 

element of the FFUX 
RADIUS vector. 


7.5.2 Kron Radius 

The Kron radius is the typical size of the aperture already described in the AUTO photometry 
section (section l7.4.3b . The whole photometry process is controlled by the PHOT_AUTOPARAMS 
and the definition of the Kron radius is given in section 1X5.II 


Parameter 

Default 

Type 

Description 

PHOT_AUTOPARAMS 

2.5, 1.5 

floats (n=2) 

MAG_AUTO 

parameters: 

(Kron_fact),(min_radius) 

PHOT_AUTOAPERS 

0.0,0.0 

floats (n = 
2) 

MAG AUTO min¬ 
imum (circular) 

aperture diameters: 
estimation disk, and 
measurement disk. 


The Kron_fact is the numbe of Kron radii the aperture is set at and the min_rad is the minimum Kron 
radius for which this is done. Otherwise the minimum aperture specified in PHOT_AUTOAPERS 
is used for the photometry. 

7.5.3 Petrosian Radius 

The Petrosian radius is another defined radius for photometry and this parameter has also only one 
input parameter associated with it, PHOT_PETROPARAMS. There is however a defining param¬ 
eter for the Petrosian parameter and this is the ratio r/. Mostly this is set to 0.2 but occasionally 
it is set to 0.5. Unfortunately, it is impossible to change this at this time. See for the complete 
expression of the Petrosian radius section lX5.2l 

Parameter Default Type Description 

PHOT_PETROPARAMS 2.0,3.5 floats (n=2) MAG_PETRO pa^ 

rameters: (Pet- 

rosian_fact),(min_radius) 

The Petrosian_fact is the number of Petrosian radii the aperture is set at and the min_rad is the 
minimum Petrosian radius for which this is done. Otherwise the minimum radius is used for the 

'^Leisurely background reading on the subject does not excist of course but there is ? who proposed the radius for 
accurate galaxy photometry and more recently (by about 25 years) lOraham and Driv^ l200.5ll who gives a very nice 
overview of the whole photometry within a radius relation. Get it from astro-ph where it’s free though. You’re probably 
a Graduate student and poor. 
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photometry. 

Figure 7.7: In the checkimage there are two radii visible. Good to know: as soon as you specify 
either AUTO (Kron) parameters or PETRO parameters for the output catalog, these apertures will 
be drawn in the checkimage (APERTURES). The Petrosian is usually the outer radius. 


NOTE: The output parameters asked for and the checkimage with the apertures 
(checkimage_type = APERTURES) influence each other. If no parameter de¬ 
pending on the Kron radius (all the AUTO photometry and KRON_RADIUS) 
is asked for in the parameter file, the Kron radius is NOT drawn in the ckeckim- 
age. 


*’lf you are having deja-vu all over again, it’s because 
is how the SE code was made as well. 


copies this directly from the section above...which I’m sure 
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7.5.4 Masking Overlapping Objects 

Now what if there are two objects overlapping each other? How to account for the overlapping pix¬ 
els? This is handled by the MASK_TYPE parameter. NONE means that the counts in the overlap 
are simply added to the objects total. BEANK sets the overlapping pixels to zero. CORRECT, the 
default, replaces them with their counterparts symmetric to the objects’ center. Best if you leave it 
at default. I’m just mentioning it out of completeness 


Parameter Default Type 

Description 

MASK_TYPE CORRECT keyword 

Method of masking of 
neighbors for photom¬ 
etry: 

NONE 

no masking. 

BEANK 

put detected pixels be¬ 
longing to neighbors to 


zero. 

CORRECT 

replace by values of 
pixels symmetric with 
respect to the source 
center. 


7.6 SE Running 

These inputparameter govern the way SE runs, if it should heed flags, how it should heed those, if 
and what to put in an outputimage, how much it should comment and how much memory it should 
use. 

7.6.1 Flags 


Figure 7.8: Well sometimes flags are important. Cartoon by S. Bateman, used without any per¬ 
mission whatsoever but I don’t make money off this anyway. 

If pixels in your image should be flagged as unreliable or other, SE can use a flag image for this 
purpose. This is well described (I think) in the official manual so I copies that section in section 
??. The internal flags of SE are described in ??. If however you have some kind of quality image 
(such as a Drizzle weight image or a coverage map coming out of the MOPEX pipeline), then you 
could conceivably convert this to a flag image to be fed to SE here. SE will then combine your 
flags (from the flagimage^^) with it’s own internal flags^°. It also has different ways to combine the 

'*N0NE might be useful to get the total flux from a large extended underlying structure with bright patches and no 
foreground stars. 

‘■^Specified in FLAGJMAGE...ta-dah! 

^°The ones that state that Timmy...I mean your object is too close to the edge etc etc... 
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flags (specified in FLAG_TYPE). If your flags are in ascending order of awfulness (flag = 1 means 
okay but with a bad pixel, flag=100 means you made this part of the image up...) then you could 
go for MAX or OR option 


Parameter 

Default 

Type 

Description 

FITS.UNSIGNED 

N 

boolean 

Force 16-bit FITS input data to 
be interpreted as unsigned inte¬ 
gers. 

FEAGJMAGE 

flag.fits 

strings 
(n < 4) 

File name(s) of the flagimage(s) . 

FEAG.TYPE 

OR 

keyword 

Combination method for flags on 
the same object: 



OR 

arithmetical OR, 



AND 

arithmetical AND, 



MIN 

minimum of all flag values. 



MAX 

maximum of all flag values. 



MOST 

most common flag value. 


7.6.2 Interpolation 


If the data for pixels is missing, SE can interpolate. These parameters regulate the interpolation. 
Best kept at default. Eor some, the x and y gaps allowed are a bit wide (16 pixels after all, it is 
almost an entire object...). On the other hand, it allowd Sextractor to give you catalogs despite bad 
columns. Therefore do not set to zero. 


Parameter 


Default Type 


Description 


INTERP_MAXXFAG 

16 

integers 

Maximum x gap (in pixels) al¬ 



(n < 2) 

lowed in interpolating the input 
image(s). 

INTERPJMAXYFAG 

16 

integers 

Maximum y gap (in pixels) al¬ 



(n < 2) 

lowed in interpolating the input 
image(s). 

INTERP.TYPE 

AFF 

keywords 

Interpolation method from the 



(n < 2) 

variance-map(s) (or weight- 
map(s)): 



NONE 

no interpolation. 



VAR_ONFY interpolate only the variance- 




map (detection threshold). 



AFF 

interpolate both the variance- 
map and the image itself. 
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7.6.3 Memory Use 

These are the parameters regulating the memory use of SE. To be honest, they are best kept at the 
values in the ./sextractor2.2.2/config/default.sex file. If you have the FEAGS output, then you 
ean eheek if there were any memory problems in the SE run. If so, this ean be very bad for the 
eompleteness of your eatalog (it isn’t probably). In that ease you might want to fiddle with the 
MEMORY_BUFSIZE and rerun. 


Parameter 

Default Type 

Deseription 

MEMORY_BUFSIZE 

integer 

Number of seanlines in the im- 
agebuffer. Multiply by 4 the 
frame width to get equivalent 
memory spaee in bytes. 

MEMORY.OBJSTACK 

integer 

Maximum number of objeets 
that the objeetstaek ean eon- 
tain. Multiply by 300 to get 
equivalent memory spaee in 
bytes. 

MEMORY_PIXSTACK 

integer 

Maximum number of pixels 
that the pixel-staek ean eon- 
tain. Multiply by 16 to 32 to 
get equivalent memory spaee in 
bytes. 


7.6.4 Neural Network 

There is to date only one neural network file and it’s in the same direetory (default.nnw) as the other 
eonfig files. Use this one. Don’t edit it. I think the original idea was to have speeialized neural 
network files for different types of instruments but it turns out it’s mueh easier to run something on 
a SE eatalog. 

Parameter Default Type Deseription 

STARNNW_NAME - string Name of the file eontaining 

the neural network weights for 
star/galaxy separation. 


7.6.5 Comments 

The VERBOSE_TYPE parameter regulates the amount of eomments printed on the eommand line. 
It eould possibly be instruetive to run it with FUEE onee in a while, the deseriptions are not very 
helpful but then again you only want to use QUIET if SE is part of some kind of pipeline and FUEE 
if something is off and you ean’t figure out what. 

^*SE was programmed to do large images, even with limited memory and computing power. The only reason to 
change these defaults is when you get something of a stack overflow. 
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Parameter Default Type Deseription 

VERBOSE.TYPE NORMAL 

keyword How mueh SExtraetor com¬ 

ments its operations: 

QUIET run silently, 

NORMAL display warnings and limited 

info concerning the work in 
progress, 

EXTRA_WARNINGS 

like NORMAL, plus a few 
more warnings if necessary, 
LULL display a more complete infor¬ 

mation and the principal pa¬ 
rameters of all the objects ex¬ 
tracted. 


7.7 SE output settings 

SE has two types of output. The catalogs with a whole range of characteristics of each of the 
detected objects and outputs which allow you to compare SE estimates of background, apertures 
and objects with the real data. 

7.7.1 Catalog 

The catalog is what you are running SE for! So in CATALOG_NAME , you specify the name of 
the output catalog. Again it’s probably a good idea to start straight away with a naming convention, 
(a .cat extension for instance!) 

The CATALOG_TYPE enables you to specify the type of outputcatalog. Personally I prefer the 
ASCIIJ4EAD, as it allows me to read it in just about anywhere and still tells me which parameters 
are listed. The nice thing about the fits catalog is that all the input parameter settings are saved in 
the header. 


NOTE: the fits option can’t handle array outputinformation such as 
MAG_APER,ELUX_RADIUS if more than one value! 


22 

The ASCILSKYCAT option for instance does not list all the parameters before the actual 
catalog (like ASCIIJ4EAD) but puts the name of the output parameter on top of the column in 
question. 

And which parameters to list is specified in fhe file given fo PARAMETERS_NAME. 

have this on Ed’s authority. I have NO idea how to display a fits table. Since there have been quite a number of 
updates, I think this may be fixed by now. Not very helpful I know. 
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Parameter Default Type Deseription 

CATALOG_NAME - string Name of the output eatalog. If 

the name “STDOUT” is given 
and CATAEOG TYPE is set to 
ASCII, ASCII HEAD, or ASCII 
SKYCAT, the eatalog will be 
piped to the standard output (std- 
out) 

CATAEOG_TYPE - keyword Format of output eatalog: 

ASCII ASCII table; the simplest, but 

spaee and time eonsuming, 

ASCII_HEAD 

as ASCII, preeeded by a header 
eontaining information about the 
eontent, 

ASCILSKYCAT 

SkyCat ASCII format (WCS eo- 
ordinates required), 

FITS_1.0 FITS format as in SExtraetor 1, 

FITSXDAC 

FITS “EDAC” format (the origi¬ 
nal image header is eopied). 

PARAMETERS _NAME - string The name of the file eontaining 

the list of parameters that will be 
eomputed and put in the eatalog 
for eaeh objeet. 


7.7.2 ASSOC parameters 

These are the parameters dealing with erosseorrolating two eatalogs: a eatalog of targets and the 
output eatalog. (The target catalog is given in ASSOC_NAME, the output one is created as SE 
runs). The cross correlation is controlled by two parameter: ASSOC_RADIUS and ASSOC_TYPE, 
the first governing the search radius and latter which objects gets selected if there are multiple 
candidates near the positions. The numbers of the columns in target catalog which contain the 
x,y positions of the objects need to be in ASSOC_PARAMS. A third column here can be used 
as weight. So you can crosscorellate catalogs weighted with magnitude but also with another 
parameter such as FWHM or no of pixels. Useful it you’re crosscorellating between catalogs of 
different filters. Some of the columns in target catalog (with the to-be-crosscorellated objects) can 
be put into the second output catalog. These columns are specified in ASSOCJDATA and end up 
in fhe VECTOR_ASSOC. 




SE OUTPUT SETTINGS 


45 


NOTE: ASSOC works only with pixel positions (NOT RA and DEC!) 


NOTE: Be aware of shifts and rotation between images when erosseorrellating 
two eatalogs. Eind shifts beween images with imeentroid in IRAE for instanee. 


NOTE: ASSOC will appear not to work if you don’t ask for a ASSOC output 
parameter in the parameter file. Your eatalog better eontain either the VEC- 
TOR_ASSOC or the NUMBER_ASSOC output, otherwise SE simply runs and 
outputs all objeets it has deteeted, not just the erosseorellated ones^^. 


Emmanuel sez: Unfortunately, SExtraetor ean only handle floating point num¬ 
bers in ASSOC files. In fael, ASSOC file may have eommenf lines (like an SE 
header) buf no lab spaeing eilher. Your besl bel is jusl a lisl of x and y positions 
(wilh possibly a weighl) separafed wifh while spaees and nolhing else. 
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Parameter 

Default 

Type 

Deseription 

ASSOC_NAME 

sky.list 

string 

Name of the ASSOC 

ASCII file. 

ASSOC_PARAMS 

2,3,4 

integers (n < 

Nos of the columns 



2, n < 3) 

in the ASSOC file 
that will be used 

as coordinates and 
weight for cross¬ 
matching. 

ASSOC_RADIUS 

2.0 

float 

Search radius (in pix¬ 
els) for ASSOC. 

ASSOC.TYPE 

MAG_SUM 

keyword 

Method for cross¬ 
matching in ASSOC: 



FIRST 

keep values corre¬ 
sponding to the first 
match found. 



NEAREST 

values corresponding 
to the nearest match 
found. 



MEAN 

weighted-average 

values. 



MAGJMEAN 

exponentially 
weighted average 

values. 



SUM 

sum values. 



MAG_SUM 

exponentially sum 
values. 



MIN 

keep values cor¬ 
responding to the 
match with minimum 
weight, 



MAX 

keep values cor¬ 
responding to the 
match with maxi¬ 
mum weight. 
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Paramefer 

Defaulf 

Type 

Descripfion 

ASSOCSEEEC.TYPE 

MATCHED 

keyword 

Whaf sources are 
prinfed in fhe ouf-puf 
cafalog in case of 
ASSOC: 



AEE 

all defections. 



MATCHED 

only mafched defec¬ 
tions. 



-MATCHED 

only defections fhaf 
were nof mafched. 

ASSOC_DATA 

2,3,4 

infegers 

Numbers of fhe 



(n < 32) 

columns in fhe AS¬ 
SOC file fhaf will be 
copied fo fhe cafalog 
ouf-puf. 


NOTE: ASSOC does not work as long as you do not ask for the ASSOC output 
in the parameter file (specified in PARAMETERSJSAME.). So don’f forgef 
fo ask for VECTORJ^SSOC or NUMBER_ASSOC in fhe paramefer file or 
ofherwise all fhe defecfions are reporfed. 


7.7.3 The Check-images 

SE can oufpuf some of fhe maps used in infermediafe sfeps. The names for fhe oufpuf fifs files 
are specified in CHECKIMAGEJSAME. Keep a convenfion like your file_background.fifs or your 
file_bgr.fifs. A lisf of up fo 16 can be given (separafed by a comma). 

The fype of oufpuf files you wanf is defined in CHECKIMAGE_TYPE. As you can see, mosf 
of fhese have fo do wifh fhe background esfimafion. Nofable excepfions are fhe APERTURES and 
SEGMENTATION options. APERTURES is a good diagnosfic on whefher or nof your fhreshold is 
righf and fhe SEGMENTATION will show you if fhe objecfs are broken up too much or nof. Eoad 
fhe original and fhis segmenfafion image info saofng and compare. 

SEGMENTATION has anofher useful feafure, fhe number in fhe cafalog is given as fhe value 
fo fhe isoarea in fhis image. Good for figuring ouf whaf is whaf from fhe cafalogs. Also, fhe 
segmenfafion image is used as inpuf for follow-up fif programs of exfended sources (GAEEIT and 
GIM2D). 


^"^Oddly enough, this more than covers all your options, including NONE for the type of checkimage. 
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NOTE: APERTURES checkimage shows the apertures you’ve asked for. So if 
you ask for a MAG_AUTO, the Kron radius will he drawn on the eheekimage. 
If NO Kron derived parameter is put into the eatalog, the Kron radius will 
also not be drawn on the APERTURES eheekimage. Same for the APER and 
PETRO apertures. 


Parameter Default Type Deseription 

CHECKIMAGE_NAME eheek.fits strings (n < File name for eaeh 

16) eheek-image . 
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Parameter 

CHECKIMAGE.TYPE 


Type _ 

keywords (n < 16) 

NONE 

IDENTICAE 

BACKGROUND 

BACKGROUND_RMS 

MINIBACKGROUND 

MINIBACK_RMS 

-BACKGROUND 

EIETERED 

OBJECTS 

-OBJECTS 

APERTURES 

SEGMENTATION 


Deseription 
Type of information 
to put in the eheek- 
images : 
no eheek-image, 
identieal to input im¬ 
age (useful for eonvert- 
ing formats), 
full-resolution interpo¬ 
lated baekground map, 
full-resolution interpo¬ 
lated baekground noise 
map, 

low-resolution baek¬ 
ground map, 
low-resolution baek¬ 
ground noise map, 
baekground-subtraeted 
image, 

baekground-subtraeted 
filtered image (requires 
PIETER = Y), 
deteeted objeets, 
baekground-subtraeted 
image with deteeted 
objeets blanked, 

MAG APER and MAG 
AUTO integration lim¬ 
its, 

display patehes eorre- 
sponding to pixels at¬ 
tributed to eaeh objeet. 


Figure 7,9: Illustration of the different eheekimages possible. The original inputimage, the BACK¬ 
GROUND image, the EIETERED image, the OBJECTS image, the SEGMENTATION image and 
the APERTURES image, the eontrast of the BACKGROUND image has been exaggerated. 





Chapter 8 

Output Parameters 


The catalogs with output parameters is what the whole exercise is all about! ^ You can finally start 
constructing your Hertzsprung-Russel diagrams or lensing shear fields or whafever. The param¬ 
eters you wanf in your cafalogs should be listed in fhe file you gave fo PARAMETERS _NAME. 
Unless you keep using fhe same file for fhis, I really recommend using ASCII_HEAD type cafalogs 
^ The oufpuf cafalogs will have a nice header wifh a lisf of all fhe paramefers. 

Parameters in fhe SE oulpufcafalog can be divided info geomefric paramefers ,phofomefric 
paramefers, asfromefric paramefers and fiffed parameters.^ Geomefric paramefers will fell you 
whaf shape fhe objecf is in (basically how fhe lighf of fhe objecf is disfribufed over fhe pixels of 
fhaf objecf) and fhe phofomefric paramefers fell you simply how much lighf fhere is. Asfromefric 
paramefers give fhe posifion of fhe objecf in fhe image, be if in pixels or ofher coordinates.^ Eiffed 
parameters are calculated from filling for inslance a PSE fo fhe dala of fhe objecf. Mosl of Ihese 
are slill being developed and do nol work yel.^ 

There are a few fhaf do nol fall in any of Ihese convenienl calegories; fhe calalog number, fhe 
flag paramefers and fhe parameters associated wifh crosscorrellaling cafalogs. 

Name description unil 

NUMBER Running objecf number 

EEAGS Exlraclion flags 

IMAEEAGS JSO EEAG-image flags OR’ed over fhe iso. profile 

NIMAEEAGS JSO # flagged pixels entering IMAEEAGS JSO 


'Well that and the check-images perhaps. 

^By setting the parameter CATALOG-TYPE to ASCII_HEAD in the configuration file. 

^Just kidding, SE does not fit anything yet. It’s so fast for a reason. 

"*There has been some improvement here as well. The now nearly ubiquitous World Coordinate System WCS is also 
used. 

^Meaning, I’ve tried but have not gotten them to work yet...If you have, let me know. 
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8.1 Photometric Parameters 

Photometric parameters are either flux or magnitude determined by SE. However SE has five dif¬ 
ferent ways of determining these; isophotal, isophotal-corrected, automatic, best estimate and aper¬ 
ture. These are discussed in the Photometry section in the input chapter (section lT!^ . To recap: 


ISO 

ISOCOR 

AUTO 

BEST 

APER 

PETRO 

PROEIEE 


Photometry derived from the counts above the threshold 
minus the background (see also section IT.d.H . 

ISO photometry, corrected for loss as a Gaussian profile 
(see also section lT.4.21) . 

Photometry from the Kron flexible elliptical aperture. ? 
(see also section IT.4.31) 

Choice between AUTO and ISOCOR. AUTO, except 
when influence from neighbors is more than 10%. (see 
also section IT.4. 4b 

Photometry from circular, user specified 
(PHOT_APERTURES in the config file), apertures, 
(see also section lT.4.5b 

Photometry from the Petrosian aperture, very similar to 

the Kron aperture, (see also section ??) 

the weighted photometry using the ’filtered’ image for the 

weight. 


Figure 8.1: Illustration of the different apertures possible; ISO, ISOCOR, AUTO and APER (user 
specified in PHOT_APERTURES) 

There are two other photometric parameters of interest: MUJVIAX, the surface brightness of 
the brightest pixel and the MU_THRESHOED, the the surface brightness corresponding to the 
threshold. This last parameter is good to inspect if the threshold is set with respect to the back¬ 
ground RMS. The value of the Background at the position of the object is often also interesting to 
know. Especially when deciding whether or not to switch between GEOBAE and EOCAE in the 
BACK_TYPE parameter (see section IT^ . 

The Petrosian aperture is a recent addition to SE (since v2.4.4 as far as I know) and it is very 
similat to the Kron radius. The apreture has a different radius (a Petrosain vs. a Kron one.). 
However the position angle and the ellipticity are the same as the Kron aperture. The Petrosian 
radius is usually bigger than the Kron one. See also section lT.5.2ll8.5.l1 and l7.4.3] for Kron aperture 
stuff and see section lT.5.3l and l8.5.2l for more on the Petrosian radius. 


NOTE: oh and I’ll say it again...the ’’BEST” photometry is a misnomer. Rather 
use something that is consistent across your image. 
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NOTE: For a color (Great Brittain and Canada: colour) measurement, the ISO 
and APER options are good, especially when run in dual mode since you’ll 
know the apertures are the same. The other apertures may be too inclusive in 
crowded fields. 


Name 

description 

unit 

FEUXJSO 

Isophotal flux 

count 

FEUXERRJSO 

RMS error for isophotal flux 

count 

MAGJSO 

Isophotal magnitude 

mag 

MAGERRJSO 

RMS error for isophotal magnitude 

mag 

FEUXJSOCOR 

Corrected isophotal flux 

count 

FEUXERRJSOCOR 

RMS error for corrected isophotal flux 

count 

MAGJSOCOR 

Corrected isophotal magnitude 

mag 

MAGERRJSOCOR 

RMS error for corrected isophotal magnitude 

mag 

FEUX_AUTO 

Flux within a Kron-like elliptical aperture 

count 

FEUXERR_AUTO 

RMS error for AUTO flux 

count 

MAG_AUTO 

Kron-like elliptical aperture magnitude 

mag 

MAGERR_AUTO 

RMS error for AUTO magnitude 

mag 

FEUX_BEST 

Best of FEUX_AUTO and FEUXJSOCOR 

count 

FEUXERR_BEST 

RMS error for BEST flux 

count 

MAG_BEST 

Best of MAG_AUTO and MAGJSOCOR 

mag 

MAGERR_BEST 

RMS error for MAG_BEST 

mag 

FEUX_APER 

Flux vector within fixed circular aperture(s) 

count 

FEUXERR_APER 

RMS error vector for aperture flux(es) 

count 

MAG_APER 

Fixed aperture magnitude vector 

mag 

MAGERR_APER 

RMS error vector for fixed aperture mag. 

mag 

FEUX_PETRO 

Flux within a Petrosian-like elliptical aper¬ 
ture 

count 

FEUXERR_PETRO 

RMS error for PETROsian flux 

count 

MAG_PETRO 

Petrosian-like elliptical aperture magnitude 

mag 

MAGERR_PETRO 

RMS error for PETROsian magnitude 

mag 
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8.2 Profile 

The profile opfion for phofomefry is nol a much used one. The idea is fo weigh fhe fluxes wifh 
fhe values in fhe filfered (smoofhed) image (see section 17.3.21 for fhe opfions here), in general, 
smoofhing correlafes fhe noise in an image and broadens profiles. By weighting wifh fhe smoofhed 
profile, you give exfra weighl fo fhe brighfesf pixels (as opposed fo counting all fhe flux equally). 


NOTE: fhe PROFILE opfion seems fo be a fix for when you fhink you include 
loo much noise in fhe flux measuremenfs (i.e. fhe defection Ihreshold is loo 
low, background loo funky or somelhing...) Use wifh some caulion. 


And remember kids: fhe elliplicify and everylhing is derived from fhe ’tillered’ 
(i.e. smoofhed) image. The PROFILE option in Ihis conlexf mighf be a good 
way fo check up how much differenl Ihings are belween fhe smoofhed and fhe 
original image. 


Name 

description 


unif 

FLUX_PROFILE Flux weighfed by fhe FILTERed profile 

FLUXERR_PROFILE RMS error for PROFILE flux 
MAG_PROFILE Magnilude weighfed by fhe FILTERed pro 

tile 

MAGERR_PROFILE RMS error for MAG_PROFILE 

counl 

counl 

- mag 

mag 

Name 

description 

unif 


MU.THRESHOLD 

Defecfion Ihreshold above back¬ 
ground 

mag X 

arcsec—2 

MU_MAX 

Peak surface brighfness above back¬ 
ground 

mag X 

arcsec“^ 

BACKGROUND 

Background al cenlroid position 

counls 


THRESHOLD 

Defecfion Ihreshold above back¬ 
ground 

counfs 



8.3 Aslromefric Paramelers 

The aslromefric paramelers are simple enough; Ihey fell you where fhe objecl is localed. However 
which pixel do you lake for fhe cenfer of fhe objecfs? The maximum flux pixel? Or fhe barycenler? 
^ Do you wanf if in RA and DEC or in prosaic pixels? All of fhese you should be able fo specify 
from fhe lisl below. 


®The Barycenter is the flux-weighted average position, the first moment, SE computes all the second moments as 
well, see the geometric parameters section and chapter 9 in the manual v2.1.3. 
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NOTE: the x and y positions (with the exeeption of ’peak’ values) are the 
baryeenters of objeets; the weighted mean position. 


The definition for the baryeenter of an objeet is: 


X = x = 




Y = y = 


^liUi 


This is the first order moment of the objeet.^ The minima and maximum x and y pixelvalues almost 
speak for themselves. These are the most extreme values for x and y that are still within the objeet. 
Handy if you want to make postage stamps of some of your objeets from the fits file. The PEAK 
values are fhe x,y and ra and dee values for fhe posifion of fhe brighfesf pixel. The value for fhe flux 
and surfaee brighfness of fhaf pixel ean be found in oulpuf-paramefers EEUXJVIAX and MU_MAX 
respeefively. The differenee befween fhe position of fhe brighfesf pixel and fhe (bary) eenfer of 
an objeef mighf help you idenlify blended objeefs. If fhey differ mueh, fhan you’ve gol somefhing 
lobsided. 


NOTE: all fhe Xi and Ij values in fhese formulae are fhe values from fhe pix¬ 
els identified in fhe segmenfafion map as belonging fo fhe objeef. If ean be 
influeneed by defeefion fhresholds and segmenfafion settings. 


NOTE: The origin for XJMAGE efe is pixel 1,1. SO The firsl pixel in fhe lefl 
bottom eorner of fhe EITS image is number 1,1 and nof 0,0. 


^Yes this is straight from the sextractor manual. All this stuff is quite complete in the manual v2.1.3 and above. Is 
there a problem? 
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Name 

description 

unit 

XMINJMAGE 

Minimum x-coordinate among detected pix¬ 
els 

Minimum y-coordinate among detected pix¬ 
els 

Maximum x-coordinate among detected pix¬ 
els 

Maximum y-coordinate among detected pix¬ 
els 

y-coordinate of the brightest pixel 

pixel 

YMINJMAGE 

pixel 

XMAXJMAGE 

pixel 

YMAXJMAGE 

pixel 

YPEAKJMAGE 

pixel 

XPEAK_WORLD 

World-x coordinate of the brightest pixel 

deg 

YPEAK_WORLD 

World-y coordinate of the brightest pixel 

deg 

ALPHAPEAK_SKY 

Right ascension of brightest pix (native) 

deg 

DELTAPEAK^KY 

Declination of brightest pix (native) 

deg 

ALPHAPEAKJ2000 

Right ascension of brightest pix (J2000) 

deg 

DELTAPEAKJ2000 

Declination of brightest pix (J2000) 

deg 

ALPHAPEAK_B1950 

Right ascension of brightest pix (B1950) 

deg 

DELTAPEAK_B1950 

Declination of brightest pix (B1950) 

deg 

XJMAGE 

Object position along x 

pixel 

YJMAGE 

Object position along y 

pixel 

XJMAGE_DBL 

Object position along x (double precision) 

pixel 

YJMAGE_DBL 

Object position along y (double precision) 

pixel 

X.WORLD 

Barycenter position along world x axis 

deg 

Y.WORLD 

Barycenter position along world y axis 

deg 

XJMAMA 

Barycenter position along MAMA x axis 

m-6 

Y_MAMA 

Barycenter position along MAMA y axis 

m-6 

ALPHA_SKY 

Right ascension of barycenter (native) 

deg 

DELTA_SKY 

Declination of barycenter (native) 

deg 

ALPHAJ2000 

Right ascension of barycenter (J2000) 

deg 

DELTAJ2000 

Declination of barycenter (J2000) 

deg 

ALPHA_B1950 

Right ascension of barycenter (B1950) 

deg 

DELTA_B1950 

Declination of barycenter (B1950) 

deg 


According to the v2.3 manual of SE, FITS header-information is used for the WORLD coor¬ 
dinates. World Coordinates System (WCS) is often what you use to line up different images from 
different observatories in ds9 (that’s what I do...). 

No idea what the MAMA axis is. 

8.4 Geometric Parameters 

Geometric parameters describe the shape and size of the object. SE computes the moments of 
an object and determines elliptical parameters from these. Both the moments and the elliptical 
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parameters and derivatives from these can be included in the output. These parameters are treated 
much more extensive in chapter 9 of the manual (version 2.1.3). 


8.4.1 Moments 

The first order moments are the barycenters of course but the second order moments can also be 
given (see section 9.1.4 in the manual v2.1.3): 


X2 = 


^li 


— X 


( 8 . 1 ) 


Y2 = 


_2 

- - - V 

Yli ^ 


(8.2) 


XY =^ = 


Yli 


- xy 


(8.3) 


Name description unit 

X2JMAGE Variance along x pixel^ 

Y2 JMAGE Variance of position along y pixel^ 

XYJMAGE Covariance of position between X and y pixel^ 

X2_WORED Variance of position along X-WORED (alpha) deg^ 

Y2_WORED Variance of position along Y-WORED (delta) deg^ 


XY.WORED Covariance of position X-WOREDA"-WORED deg^ 


NOTE: By adding ERR+parameter, the error can be obtained as well. (e.g. 
ERRX2 JMAGE) 


8.4.2 Ellipse parameters 

Erom these moments, the position angle, the minor and major axis of an ellipse can be derives 
(section 9.1.5 in the manual v2.3) and also a second way of describing this ellipse is given (sec 
9.1.6). The official manual gives the whole derivation of these parameters but since we want a 
quick result ITl cut it down some. 


NOTE: All these elliptical parameters are computed in the SMOOTHED image. 
So If you have used a smoothing kernel in the input, then the EEEIPTICITY is 
also of a smoothed object. * 
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Minor and Major Axes 

The following definitions for these computed parameters are: (see for the derivation the v2.3 
manual page 28) 


A2 


2 


+ x^ 
2 



(8.4) 


(8.5) 


These are illustrated in figure lOl Not also that the Petrosian and Kron radii (section 18.5.21 and 
18.5.1b are expressed in the same units as those used for the A and B values. 

Figure 8.2: Illustration of the Kron radius, x and y are the computed center in of the object. 


Name 

description 

unit 

AJMAGE 

Profile RMS along major axis 

pixel 

BJMAGE 

Profile RMS along minor axis 

pixel 

THETAJMAGE 

Position angle (CCW/x) counterclockwise!!!^ 

deg 

A.WORED 

Profile RMS along major axis (world units) 

deg 

B.WORED 

Profile RMS along minor axis (world units) 

deg 


Errors can be obtained using ERR+parameter in the paramfile. 

Position Angle 

This is the rotation of the major axis with respect to NAXIS1 (x-axis...) counterclockwise. The def¬ 
inition of the position angle is different (pa is with respect to North.) However the conversion can 
obviously be done and the Position Angle (PA) is given in other THETA parameters (SKY, B1950 
and J2000). This can be very useful if you’re looking for slit-positions for your spectrograph. 


Name 

description 

unit 

THETA.WORLD 

Position angle (CCW/world-x) 

deg 

THETA_SKY 

Position angle (east of north) (native) 

deg 

THETAJ2000 

Position angle (east of north) (J2000) 

deg 

THETA_B1950 

Position angle (east of north) (B1950) 

deg 


Errors can be obtained using ERR-i-parameter in the paramfile. 

The Other Ellipse parametrisation 

So the manual also gives another parametrisation of the elliptical aperture. I have not used these 
parameters. I have no idea what to do with them. Maybe you do. In figure lOl the relation between 
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these parameters and the moments is given. 
Name description 


unit 


CXX JMAGE 

Cxx object ellipse parameter 

pixel“^ 

CYY JMAGE 

Cyy object ellipse parameter 

pixel“^ 

CXY JMAGE 

Cxy object ellipse parameter 

pixel“^ 

CXX.WORLD 

Cxx object ellipse parameter (WORLD units) 

deg-2 

CYY.WORLD 

Cyy object ellipse parameter (WORLD units) 

deg-2 

CXY.WORLD 

Cxy object ellipse parameter (WORLD units) 

deg“2 


Errors can be obtained using ERR+parameter in the paramfile. 


Ellipticity and Elongation 

How stretched the object is can be parameterised with weither the Ellipticity or the Elongation 
parameters, both of which are simply computed from the minor and major axes. 


A 

ELONGATION = - 
B 


( 8 . 6 ) 


TD 

ELLIPTICITY = 1 - - = 1 
A 


I 

ELONGATION 


(8.7) 


Name description unit 

ELONGATION A_IMAGE/B JMAGE 
ELLIPTICITY I - B_IMAGE/AJMAGE 


NOTE: If you have EILTERED the image for detection then your ellipticity and 
elongation values will be more round than the actual object. 


Only the ELONGATION and the ELLIPTICITY parameters do not come with a error. Note 
that there is an extensive discussion of these in section 9.1.8 of the manual. 

8.4.3 Area Parameters 

This can be done by the parameters describing the isophotes as fitted by SE or the conclusion of 
the neural network classification. To get an idea of the size of an object, reasonably independent of 
the brightness, the EWHM or the ELUX_RADIUS 

SE divides an object up into 7 isophotes above the ANALYSIS_THRESH. The areas above the 
isophotes is fed to the neural network. These can be put into the catalog with the ISOn parameters. 

'°rd be saying the exact same thing and typing in that many functions is a pain, even in LaTeX. 





60 


CHAPTER 8: Output Parameters 


Name 

description 

unit 

ISOAREA.WORED 

Isophotal area above Analysis threshold 

deg^ 

ISOAREAE.WORED 

Isophotal area (filtered) above Detection 
threshold 

deg2 

ISOO 

Isophotal area at level 0 

pixel^ 

ISOl 

Isophotal area at level 1 

pixel^ 

IS02 

Isophotal area at level 2 

pixel^ 

IS03 

Isophotal area at level 3 

pixel^ 

IS04 

Isophotal area at level 4 

pixel^ 

IS05 

Isophotal area at level 5 

pixel^ 

IS06 

Isophotal area at level 6 

pixel^ 

IS07 

Isophotal area at level 7 

pixel^ 


8.4.4 Full-Width Half Max 

Assuming a Gaussian profile for the object, a Full-Width at Half Maximum can be computed. No 
object is a Gaussian -especially in the wings- but this does give you a reasonable idea of the PSF 
width in the case of a star for instance. It’s a useful thing to have. 


Name description unit 


EWHM.WORED 

FWHM assuming a Gaussian core 

pixel 

EWHM.WORED 

FWHM assuming a Gaussian core 

deg 

VIGNET 

Pixel data around detection 

count 

VIGNET_SHIET 

Pixel data around detection corrected for 

shift 

count 

THRESHOEDMAX 

Maximum threshold possible for detection 

count 


8.5 Radii 

Indicators for the size of objects are multiple and SE can provide a lot of them if so desired. 
Especially for extended source several schemes have been developed to define an aperture to deter¬ 
mine the total flux and characterize the size of an object. In the section 17.5.2118.4.21 and 17.4.^ the 
Kron aperture and the various associated parameters were already discussed. The KRON_RADIUS 
output parameter is the indicator of size of the KRON aperture. Alternatively, there is the PET¬ 
ROSIAN.RADIUS (since SE v2.4.4). Both of these are expressed in multiples of major axis B. 

8.5.1 Kron radius 

? defined the KRON radius to get 90% of an objects light as follows: 

2'k I{x)x^dx 
2t: ff^I(x)xdx 




( 8 . 8 ) 












Radii 


61 


Table 8.1: Radii nomenclature 


Radius 

symbol 

Description 

Kron 

Ri 

The typical size of the flexible aperture -computed from moments- 
defined by ?. 

Petrosain 

Rp 

The radius at which the surface brightness of the isophote is 

T] times the average surfacebrightness within this isophote. 

Effective 

Re 

The radius containing 50% of the total flux of an object. 

Half-light 


same radius. 

Typical 

Rt or h 

The scale in the exponential disk I = /q exp{—R/Rt) 

deVaucouleur 

R 25 

The radius at which the B-band isophote is 

25 mag arcsec“^. This is a common Radius as 
it is listed in the third reference catalog (RC3) 


According iGraham and Driven (l2005h . the orginal SE publication has changed the definition a 
bit. iBertin and AmoutsI (Il996il d^ne the Kron radius as: 


Ri 


j:ri{r) 


(8.9) 


over the two dimentional aperture and not for a lightprofile (which was Kron’s original idea.) In 
the case of a lightprofile, the Kron-radius should be according to Kron’s original definition: 


I{R) 
Si? I{R) 


( 8 . 10 ) 


So the Kron radius coming out of SE is maybe not a proper one but we’re using this as the 
estimate anyway. And if you are interested in galaxy profiles, here is a neat trick to relate this Kron 
radius to the ejfective radius (Re), the radius that encloses half the object’s light. 


Ri{x, n) = 


Re 7(3n,x) 


( 8 . 11 ) 


7 ( 2 n, x) 

In which x = b{RlRp)^/'^ and 7 is th e incomplete gamma function^^. This is all too much of a 
headache so iGraham and Driven (l2005n also give a nice table converting these radii depending on 
Sersic profile index (copied in table l8.2l) . 

Name description unit 


KRON RADIUS Kron radius in units of A or B 


And remember that the KRON flux measurements can be found in the AUTO output-parameters. 


13 


"ah these radii to characterize galaxy profiles can be a tad confusing. See table lS.ll for the definitions of all. 

*^Oh Gamma function, is there something you’re not useful for? 

'^Named differently for your inconvenience. Please remember that a lot of SE stuff got tacked on later which explains 
the funky naming ’convention’. 
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Table 8.2: Theoretical Kron Radii and Magnitudes from iGraham and Driver (2005) 
Sersic n Ri L{< 2Ri) L{< 2.5Ri) 


(Re) % % 


0.5 

1.06 

95.7 

99.3 

1.0 

E19 

90.8 

96.0 

2.0 

1.48 

87.5 

92.2 

3.0 

1.84 

86.9 

90.8 

4.0 

2.29 

87.0 

90.4 

5.0 

2.84 

87.5 

90.5 

6.0 

3.53 

88.1 

90.7 

7.0 

4.38 

88.7 

9E0 

8.0 

5.44 

89.3 

9E4 

9.0 

6.76 

90.0 

9E9 

10.0 

8.39 

90.6 

92.3 


8.5.2 Petrosian radius 


PetrosianI (119761) defined a point in the radial light profile at which the isophote at that radius was a 
certain fraction of the average surface brightness within that radius. 


2^ j^I{R')R'dR' _ {I)j, 


r]{R) = 


( 8 . 12 ) 


7rR‘^I{R) I{R) 

where the parameter rj is the fraction. This is often either 0.2 or 0.5 with 0.2 the most commonly 
used. I am pretty sure the 0.2 value is used in SE. 


NOTE: The rj cannot be changed using some input-parameter. Bummer. 


There is an input parameter associated with this in source extractor ( PHOT_PETROPARAMS) 
where the first input defines the number of Petrosian radii for the magnitude and the second the 
minimum radius. The following output parameters are associated with the Petrosian radius: 

Name description unit 

PETRO-RADIUS Petrosian apertures in units of A or B 


PEUX_PETRO 

Elux within Petrosian radius elliptical aper¬ 
ture 

count 

PEUX_PETRO 

Elux within Petrosian radius elliptical aper¬ 
ture 

count 

EEUXERR-PETRO 

RMS error for PETRO flux 

count 

MAG_PETRO 

magnitude within the N x Rp 

mag 

MAGERR_PETRO 

RMS error for the MAG_PETRO magnitude 

mag 
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More on the Petrosian radius can 


be found in iGraham and Dri^ ( 2005 ) and Graham et al. ( 2005 ). 


8.5.3 Effective radius 

The effective radius is the term commonly used to define the point in a light profile within which 
encloses half the flux from an object. In the case of SE, several radii can be defined (input parameter 
PHOT_EEUXERAC). PHOT_PEUXPRAC 0.5 is the effective radius (or SE’s determination of it.) 
Name description unit 

PEUX_RADIUS Radius enclosing a specified fraction of the pixel, 
flux 

With the photometry in the APER parameters. 

So here is another trick. If you’s want an estimate of the Sersic profile index n, you could use 
the ratio between the EEUX_RADIUS with PHOT_EEUXERAC set to 0.5 and the KRON_RADIUS 
in your catalog. Of course you would need to run some follow-up program to properly determine 
the actual index from a fit, but that is what SE is all about. 


8.6 Object classification 

This is a section completely devoted to the CEASS_STAR parameter; SE’s classification of the 
objects on the basis of a Neural Network Output. 

Name description unit 

CEASS STAR S/G classifier output none 

It can have a value between 0 (galaxy, more to the point: non-star) and 1 (star). 

8.6.1 Input Dependency 

Eortunately, you do not need to understand Neural Networks to use this but there are several input 
parameters which are directly linked to the CEASS_STAR parameter: 

• PIXEE_SCAEE: Pixel size in arcsec. (for surface brightness parameters, EWHM and star/ 
galaxy separation only). 

• SEEING_PWHM: EWHM of stellar images in arcsec. This quantity is used only for the 
neural network star/galaxy separation as expressed in the 

CEASS_STAR output. 

These are obvious. The NNW has to take in account the scale and blurring of the objects before 
judgment. These are inputparameters of the NNW. The ratio between the two is something you can 

'"^A quick and possibly dirty first look, quick sample selection etc... 

*^If you DO want to know more on Neural Networks and how they are used for complex issues such as object 
classifiacation, I found the book ’An introduction to Neural Networks’ by Kevin Gurney very useful. 
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play around with, (remember that SE does not use the pixelsseale for anything else...) But there 
are parameters where it depend more indireetly on: 

• BACK_SIZE Size, or Width, Height (in pixels) of a baekground mesh. 

• THRESH.TYPE Meaning of the DETECT THRESH and ANAEYS IS .THRESH parameters 

- REEATIVE sealing faetor to the baekground RMS. 

- ABSOEUTE absolute level (in ADUs or in surfaee brightness). 

• ANAEYSIS.THRESH Threshold (in surfaee brightness) at whieh CEASS STAR and EWHM 
operate. 

- 1 argument: relative to Baekground RMS. 

- 2 arguments: mu (magarcsec“^ ), Zero-point (mag). 


Obviously the brightness level from whieh objeets are eonsidered inliuenee heavily the elassi- 
fieation. If only the top levels are eonsidered, even elear galaxies might be elassified as stars (the 
disk lies below the threshold) and if the threshold is too low, random noise ean pass itself off as a 
faint galaxy. 


NOTE: if the threshold is relative to the RMS, the BACK_SIZE is VERY im¬ 
portant for the CEASS .STAR 


• DEBEEND JVIINCONT Minimum eontrast parameter for de-blending. 

• DEBEEND J4THRESH Number of deblending sub-thresholds. 

NOTE! if deblending is too eourse, a elump of stars ean beeome a ’galaxy’. 


NOTE! And if deblending is too pieky, a single galaxy might be ehopped up 
into several objeets. 


There is no single remedy but a DEBEEND JSTHRESH of 32 and a DEBEEND JVIINCONT 
of order 0.01 are a good plaee to start. 


HMMMM: If the ANAEYSIS.THRESH is different from the DE- 
TECT.THRESH, the objeets are NOT deteeted again. So earefull with making 
these different. 
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8.6.2 Reliability 

This section is one big note of caution. The NNW classification by SE is NOT perfect and will 
break down at the lower magnitude end. Unfortunately it does not give another value for unreliable 
classification (42 for instance) but assigns a random value between 0 and 1. (see figure lOl 

Figure 8.3: The dependence of CLASS_STAR on the luminosity of objects. The reliability clearly 
disappears at the lower end. 
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8.7 ASSOC output 

These are the two output parameters associated with cross correlating catalogs. If you have speci¬ 
fied a ASSOC file and fhe manner of cross-idenfificalion (see fhe inpuf parameters in section lOl . 
Be sure fo specify one of fhese as ofherwise fhere will be NO crosscorellafion and fhe SE cafalog 
will confain jusf as many enfries as wifhouf fhe ASSOC parameters specified.The column in fhe 
ASSOC-NAME cafalog, specified in ASSOCJDATA are given in fhe VECTOR_ASSOC 
Name descripfion unif 

VECTOR_ASSOC ASSOCiafed parameter vector 
NUMBER_ASSOC Number of ASSOCiafed IDs 


8.8 Elags 

How to steal a country through the cunning use of flags: 

Brittain: 1 claim India for Brittain! 

India: You can’t claim us! There’s 300 million of us already living here! 

Brittain: Do you have a flag? No flag, no country. 

—adapted from Eddie Izzard. 

This section is fhe one in the official manual wifh some annofafions of mine. 

A sef of bofh infernal and external flags is accessible for each objecf. Infernal flags are produced 
by fhe various defection and measuremenf processes wifhin SExfracfor; fhey fell for insfance if an 
objecf is safurafed or has been fruncafed af the edge of the image. 

External flags come from flag-maps: fhese are images wifh the same size as the one where 
objects are detected, where integer numbers can be used to flag some pixels (for instance, bad 
or noisy pixels). These types of images are becoming quite common data-products of reduction 
pipelines of current space missions like the Hubble or Spitzer telescopes or the grand surveys such 
as SEDAN. 

Different combinations of flags can be applied wifhin fhe isophofal area fhaf defines each objecf, 
fo produce a unique value fhaf will be written to fhe cafalog. 

8.8.1 Infernal Elags 

The internal flags are always computed and fhe defaull.param file lisfs fhem. Clearly They are 
accessible fhrough fhe EEAGS cafalog paramefer, which is a short infeger. EEAGS confains, coded 
in decimal, all fhe exfracfion flags as a sum of powers of 2: 

1 The objecf has neighbours, brighf and close enough to signi canfly bias fhe MAG AUTO 
phofomefry^^, or bad pixels (more than 10% of the integrated area affected). 

**I will NOT tell anyone how long it took me to figure THAT one out... 

steal the official manual’s prose from it with reckless abandon, okay? 

**This flag can only be activated when MAG_AUTO magnitudes are requested. 
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2 The object was originally blended with another one. 

4 At least one pixel of the object is saturated (or very close to). 

8 The object is truncated (too close to an image boundary). 

16 Object’s aperture data are incomplete or corrupted. 

32 Object’s isophotal data are incomplete or corrupted^^. 

64 A memory overflow occurred during deblending. 

128 A memory over flow occurred during exttaction. 

For example, an object close to an image border may have FLAGS = 16, and perhaps FLAGS 
= 8+16+32 = 56. The flags are combined such that a unique number will result as the flag entty in 
the catalogs. 

8.8.2 External Flags 

SExtractor will look for an external flags file when IMAEEAGS JSO or NIMAEEAGSJSO are 
present in the catalog parameter file. The file is specified in EEAG JMAGE The external and 
internal flags are then combined and listed as the EEAGS output parameter in the catalog. So you 
need to make sure that a combination of flags is still unique. There is a fun little exercise for you. 

It then looks for a EITS image speci ed by the EEAG IMAGE keyword in the configuration file. 
The EITS image must contain the flag-map, in the form of a 2-dimensional array of 8, 16 or 32 bits 
integers. It must have the same size as the image used for detection. Such flag-maps can be created 
using for example the WeightWatcher software (Bertin 1997).^*^ 

The flag-map values for pixels that coincide with the isophotal area of a given detected object 
are then combined, and stored in the catalog as the long integer IMAEEAGS ISO. 5 kinds of 
combination can be selected using the EEAG TYPE configuration keyword: 

OR: the result is an arithmetic (bit-to-bit) OR of flag-map pixels. 

AND: the result is an arithmetic (bit-to-bit) AND of non-zero flag-map pixels. 

MIN: the result is the minimum of the (signed) flag-map pixels. 

MAX: the result is the maximum of the (signed) flag-map pixels. 

MOST: the result is the most frequent non-zero flag-map pixel-value. 

'^An old flag inherited from SExtractor V 1.0. It has been kept for compatibility reasons (and why tear it out?). With 
SExtractor V2.0+, having this flag activated doesn’t have any consequence for the extracted parameters. 

^°SEXtractor? Weightwatchers? What’s with the naming convention here? I’m almost certain that this contributed to 
my celery-phobia. 
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The NIMAFLAGS ISO catalog parameter contains a number of relevant flag-map pixels: the 
number of non-zero flag-map pixels in the case of an OR or AND FLAG_TYPE, or the number of 
pixels with value IMAFLAGS ISO if the FLAG TYPE is MIN, MAX or MOST. 

8.9 Eitted Parameters 

There is a experimental section of SE which will hopefully become workable soon enough. This 
deals with the fitting of the Point Spread Eunction to stars and light profiles to extended objects.In 
the meantime the follow-up program GALEIT and GIM2D have cornered this market. As a result 
the priority on finishing these parameters may have slipped some... 

8.9.1 PSE fitting 

The point spread function of an instrument describes how the light from a point source is distributed 
over the detection element of a CCD. If you have a little fits file of a model of the PSE and you’ve 
figured out how to feed that to SE (if you do let me know...) then a part of SE known as PSEeX can 
fit this to every object. The following parameters should then be available to you. As you can see 
some of these are either photometric, geometric or positional parameters. I’ve just list them here 
as this is still an experimental bit of SE but should be very useful later. 

The PSE is also used for fitting the components of a galaxy in the next section. 


Name 

description 

unit 

XPSEJMAGE 

X coordinate from PSE-fitting 

pixel 

YPSEJMAGE 

Y coordinate from PSE-fitting 

pixel 

XPSE.WORLD 

PSE position along world x axis 

deg 

YPSE.WORLD 

PSE position along world y axis 

deg 

ALPHAPSE_SKY 

Right ascension of the fitted PSE (native) 

deg 

DELTAPSE^KY 

Declination of the fitted PSE (native) 

deg 

ALPHAPSEJ2000 

Right ascension of the fitted PSE (J2000) 

deg 

DELTAPSEJ2000 

Declination of the fitted PSE (J2000) 

deg 

ALPHAPSE_B1950 

Right ascension of the fitted PSE (B1950) 

deg 

DELTAPSE_B1950 

Declination of the fitted PSE (B1950) 

deg 

ELUX_PSE 

Elux from PSE-fitting 

count 

ELUXERR_PSE 

RMS flux error for PSE-fitting 

count 

MAG_PSE 

Magnitude from PSE-fitting 

mag 

MAGERR_PSE 

RMS magnitude error from PSE-fitting 

mag 

NITER_PSE 

Number of iterations for PSE-fitting 

- 

CHI2_PSE 

Reduced chi2 from PSE-fitting 

- 


All these have ERR-i-parameter uncertainties associated with them, except NITER_PSE and CHI2_PSE 
of course. That is...if these ever work... 

The reason you’d want a decent fit of the PSE to an object is to do really good photometry on 
it as you’re now accounting for al that leaked light. Note that such a correction would only work if 
the object was originally a pointsource. This is fairly useless for extended objects but for that we 
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have the galaxy fitting parameters. 


8.9.2 Galaxy profile filling 

Anolher experimenfal bil of SE. As you mighl know, average radial lighlprofile of a galaxy can be 
described by fil. In Ihe case of a spiral galaxy if is an exponential for Ihe disk and a deVaucouleur 
(r“^) profile for Ihe bulge. The (r“4) profile fils an elliplical galaxy quite well. So if Ihe galaxy is 
reasonably resolved, Ibis new bil of SE can fil Ihese profiles. If can give you Ihe bulge disk ralio’s 
and everylhing. Provided if’d aclually work.^^ Did I menlion GAEEIT and GIM2D in Ibis conlexl? 


Name 

description 

unil 

EEUX.GAEEIT 

Elux derived from Ihe galaxy fil 

counl 

EEUXERR GAEEIT 

RMS error for GAEEIT flux 

counl 

MAG.GAEEIT 

Magnilude derived from galaxy fil 

mag 

MAGERR.GAEEIT 

Magnilude error derived from galaxy fil 

mag 

ERROR.GAEEIT 

Reduced Chi-square error of Ihe galaxy 

fil 

GAEDANGJMAGE 

Galaxy disk position angle from Ihe galaxy 
fil 

Galaxy disk-scale from Ihe galaxy fil 

deg 

GAEDSCAEEJMAGE 

GAEDASPECJMAGE 

GAEDEIJMAGE 

pixel 

Galaxy disk aspecf ratio from Ihe galaxy fil 
Galaxy disk elliplicily nrl from Ihe galaxy 
fil 

Galaxy disk elliplicily nr2 from Ihe galaxy 
fil 

Galaxy bulge ratio from Ihe galaxy fil 

GAEDE2JMAGE 


GAEBRATIOJMAGE 


GAEBANGJMAGE 

Galaxy bulge position angle from Ihe 
galaxy fil 

deg 

GAEB SC AEE JM AGE 

Galaxy bulge-scale from Ihe galaxy fil 

pixel 

GAEBASPECJMAGE 

Galaxy bulge aspecf ralio from Ihe galaxy 
fil 



8.10 Principle Componenf 

Again Ihese are relaled lo Ihe psf filling SE is -hopefully- capable of in Ihe near fulure. 


As per v2.4.4 they do not... 
^^Meaning: ’possibly in your lifetime’. 
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Name 

description 

unif 


X2PCJMAGE 

PC variance along x 

pixeP 


Y2PCJMAGE 

PC variance along y 

pixel^ 


XYPCJMAGE 

PC covariance between x and y 

pixel^ 


APCJMAGE 

PC profile RMS along major axis 

pixel 


BPCJMAGE 

PC profile RMS along minor axis 

pixel 


THETAPCJMAGE 

PC position angle (CCW/x) 

deg 


PC 

Principal componenfs 





Chapter 9 

Strategies for SE use 


Baldrick: ’’Sir, I have a cunning plan ” 

Blackadder: ”As cunning as a fox who is professor of cunning at Oxford University?” 

This section is reserved for some lecturing on what I found were good tricks to use with SE. 
There are some tricks you can do to extend the depth of your exposures, to make it easier to use SE 
on a batch of objects or tune your detections more to your needs. 

9.1 Image types to use? 

There are several options open to you for the image type. It is quite usual for astronomical images 
to be the result of stacked exposures. There are several ways to stack. The interaction between the 
type of addition and the photometric inputparameters is discussed here. 

You may want to ask yourself whether you want to use the original integrated image with the 
total counts of every object or to use these divided by the exposure time, essentially counts-per- 
second images. Using counts-per-second images is not as dumb as it sounds. By using counts- 
per-second images, the dynamic range is the same for images with wildly varying exposure times. 
When examining them, it makes them much easier to compare. The signal-to-noise ratio (S/N) is 
of course still the same. 

In the case of counts-per-second images, you can determine the MAG_ZEROPOINT for an 1 
sec exposure and put this in the parameter file. Only the GAIN varies with every different exposure. 
The value for the GAIN in this case is 


Effective Gain 

Magnitude zeropoint 

Type of image 

gain X total exposure time 

zeropoint(l sec) 

input image is c/s 

gain 

zeropoint(l sec) -i- 2.5 
logioCexp. time) 

sum of N frames 

Nxgain 

zeropoint(l sec) -i- 2.5 
logio(av. exp. time) 

average of N frames 

2 xNx gain/3 

zeropoint(l sec) -i- 2.5 
logio(av. exp. time) 

median of N frames 
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the ccd gain times the exposure time. ^ 

9.1.1 Thresholds 

This last sheme works best when your thresholds are RMS noise related. Otherwise you’d have to 
modify the thresholds anyway to account for the signal-to-noise. The zeropoints of the thresholds 
remain the same for counts-per-second images and integrated images. This can be very confusing, 
using the zeropoint-for-one-sec-exposure in the MAG_ZEROPOINT parameter and the zeropoint- 
for-one-sec-exposure -i- 2.51og(exptime) as the zeropoint for the detection/analyse threshold. 

9.2 How to get faint objects? 

Eaint objects are the ahrdest and often the most interesting objects in you image. So how to get 
them? There are tricks for the detection images, the thresholds and the filters you use. 

9.2.1 Different Detection Images 

There are several strategies employed to detect and classify sources, using Source Extractor. Source 
Extractor has the in-built ability to detect sources in one image and subsequently do photometry 
on the found apertures in another image. This has the advantage that the photometry on an object 
has the same apertures in all bands and the catalogs need not be matched (the numbering for all the 
catalogs is the same). Many schemes use this ability. 


“Meta”-images 

One idea is to add data out of two or several bands and do the detection there while doing the 
photometry in the actual data. Th e image where the detection is done is referred to as a ’’meta¬ 
image” (i.e. Williams et al. ( 1996li used an I-i-V image to detect sources in the HDE). This is quite 
widespread use; you can use a median, mean or a stacked images. Just beware of differences in 
PSE when making these images. ^ 

The advantages are that apertures can be determined more accurate for objects, espescially 
faint ones. However, strucural parameters (all those computed from the object’s moments) are 
determined for the meta-image. These may not always be very representative of your objects. 


Maximum likelyhood images 

Anot her approach using m ulti-band images is to construct a image and do the detection on 
these. Szalav et al. ( 1999ah introduced this technique on the HDE. However this method appears to 
be most successful in the truly multiband (more than three at least) exposures. 

Both of these are aimed at getting as faint as possible sources. 


*As you’ve probably guessed I used counts-per-second images. It’s personal preference but it’s much easier to run 
SE in batches if the only thing you have to change is the GAIN. The zeropoint changes as soon you change filter of use 
a different ccd in an array like changing from WF2 to WF3 of the Wide Field camera on Hubhle. 

^The PSF of the resulting image should be about that of the worst PSF component. Some convolution may be needed. 


















How TO GET GOOD COLOURS OF OBJECTS? 
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Optimal image 

Most commonly, the band where your target objects are the brightest (i.e. the reddest band avail¬ 
able for field galaxies) or the band with the longest exposure is used as the detection image and the 
apertures of that band are used for the photometry in all the other bands. The detection of back¬ 
ground galaxies is done in the band with the longest wavelength available is used on the premises 
that background galaxies appear red. 

Alternatively SE can be run on each band separately and the catalogs are subsequently matches 
by position. Fixed aperture photometry, corrected for the PSF can then be used to avoid differences 
due to different apertures. This had obvious drawbacks for extended or faint sources; you don’t 
know how well the apertures match. 

9.2.2 Thresholds and Filters 

Again the settings within sextractor can be optimized for faint objects. The threshold can be set 
as low as you want it but a DETECT_THRESH of 1 ct RMS above the background is probably 
the best you can do. However how that background is determined is very imposrtant. Depending 
on your images and objects, changing the BACK_SIZE background mesh size can improve your 
detection of faint stuff, especially near bright objects (v hard in any case). Which brings us to the 
CEEAN option you might want to turn off so you can better detect faint stuff near brighter things. 

And some of the filters used to smoothe the image before detection are influential in the faint 
sources detection. Have look at the filtering section earlier. Usually a Gauss approximating the 
PSF in size is used. 


NOTE!; whatever your settings, MAG_AUTO is the best estimate of the 
magnitude but it still underestimates by as much as several tens of magnitude 
for faint objects. ^ 


Best way to figure fhis ouf is fo injecf you images wifh simulafed objecfs of known magnifude 
and check fheir magnifude in fhe SE cafalog. For simulafed objecf, I use mkobject in IRAF. Give if 
a lisf of posifions, magnifudes and profiles and you’ve gol your calibrafion objecfs. 

9.3 How fo gef good colours of objecfs? 

Colours of objecfs found in surveys are a very imporfanf fool for asfronomers. So having a good 
sfrafegy how fo gef fhe colours from SE cafalogs is probably a good idea. There are several ways fo 
defermine fhese. Do you use a defection image and phofomefry image(s) or use values from several 
independenf cafalogs? Which magnifudes fo use;ISO, ISOCOR, MAG_AUTO, MAG_APER or 
somefhing else? 
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9.3.1 separate deteetion and photometry images? 

SE has this ability to use the apertures from one deteetion image in the other, photometry image. 
This should give aeeurate eolours as the apertures the same right? Some things to eonsider: 

• The images better be aligned to the pixel. Small shifts ean result in dramatieally different 
eolours. 

• The PSE’s and seeing better be similar (i.e. the pixels in the deteetion image are indeed the 
same part of the objeet in the photometry image) 

9.3.2 whieh output to use? 

There is only two to eonsider MAG JSO and MAG_APER, the others are too dependent on the SE 
settings in the deteetion. 

MAG JSO 

The best aperture for eolours is the MAG JSO. This gives the flux in an as big an aperture that will 
fit in the objeet, will be roughly the same shape as the objeet and when using separate deteetion and 
photometry images, these will give you the most aeeurate eolours. This of eourse does not hold of 
you erosseorellate eatalogs. 

MAG_APER 

MAG_APER seems ideal for the eolour of objeets; either used in separate deteetion and photometry 
images or from erosseorrellated eatalogs, the apertures are user speeified and eomparable over the 
bands used (you ean eorreet the radii for instanee with seeing...) BUT! there are a few things to 
eonsider with fixed apertures: 

• crowding; does your aperture overlap with another object? 

• objects smaller than the aperture will have less reliable colours 

• Still, are they aligned? Is this the same part of the object? 

• In the case of small apertures, is the colour of the center of the object the same as the for the 
whole? 

• PSE aperture corrected? 

But you can take these in account. The only (dis)advantage ^ MAG_APER has is that really the 
fact that it is fixed for every object. However it is not so dependent on the detection by sextractor. 

In conclusion: with detection and photometry images, use MAG JSO, with separate detections, 
use MAG_APER. Ponder which ones would be best for you type of objects. 

"*§ 00(1 or bad? your call. 



Finding your objects of interest 
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9.4 Finding your objects of interest 

This is a very broad query. Literate searches and some intense staring at typical objects should 
give you an idea what the characteristic values of SE output of your objects are. Sometimes a 
specialised parameter added to SE might be able to help or a ratio of SE parameters (aperture 
fluxes for instance). Best thing to do is get a trainingset and start plotting. EWHM, EEEIPTICFTY, 
MAG_APER(S) PEUX_RADIUS and concentration (MUJMAX over MAG_AUTO is a reasonable 
indication but there are more) are good places to start. 

And of course colours always help. 


9.5 Strategies to find galaxies in crowded fields 


To increase SE’s ability to find field galaxies in crowded fields, several strategies may be consid¬ 
ered_ 


Gonzalez et al.l(ll998h used B-I images of NGC 3664 to detect sources instead of the I images. 


This has the advantage that stars can be negated somewhat, leaving only objects of a severely 
differing color. The disadvantage is the increase in the noise, increasing the chances of spurious 
detections. And there actually has to be a B band exposure. 

An alternative approach is to remove the stars, using DAOPHOT. It is honed to find and model 
stars in crowded fields and there is the option to produce an image with the model subtracted. 
Again, this could increase the number of spurious detections because of noise. 

Alternatively SE can simply be run on the I images and photometry can be done on the V 
image. This additional information can be used to filter the catalogs for actual galaxies. As an 
additional filter, an visual inspection of the data can be employed. 

A technique to enhance faint extended sources is called unsharp masking. Originally developed 
for photographic plate, ccd images would be smoothed and the smoothed version would be sub¬ 
tracted from the original. This technique works fine for filamentary structures like shells and tidal 
tails but is unsuited for finding galaxies in crowded fields. (SE does however smooth the image 
slightly and then does the detection) 

Metaimages, the addition of two band, suitably weighed, as detection images is another popular 
option. Again this opens the possibility of many spurious detections. 








Chapter 10 

Available Packages 


Maybe it’s not a good idea to advertise the competing software but SE might not exactly be what 
you need. Better if you figure this early. At the moment there ^e several programs and ap proaches 
in use that atternpt to detect a nd classify objects: Sextrac tor [Berlin and Amouta dlQQhh . FOCAS 
Jarvis and Tyson ( 1979 , 1981 ), DAOPHOT Stetson! ( 1987 ) and self organising maps (SOM), 


All of these have been tried and tested in relatively uncrowded fields and are expected to break 
down to some extent in crowded fields. The notable exception is DAOPHOT which is honed on 
crowded stellar fields. However it is focused on accurate stellar photometry, not on detecting 
extended objects. 

FOCAS is widely rumored not to work, especially in crowded environments but still popular. I 
personally haven’t used it but SE is somewhat more recent so I chose SE over FOCAS. 

Self Organizing Maps (SOM) are still in an experimental phase and have proved themselves 
superior to SE in classifying galaxies by using an additional neural network. However they are not 
distributed in a package. If you really need detailed classification of objects (your thesis hangs on 
it being Sab galaxies) then you’d better contact the authors of papers on SOMs or write something 
yourself. But if you can think up a reasonable parameter that will help you classify, you can always 
define your own in SE (see 

This leaves however SE as the current program of choice, almost by default. 














Chapter 11 

Follow-up Programs 


SE is now used as either a quick look or -mostly- as the first step in a more involved analysis 
that often involves fits to the images. In the case of galaxies, there are several packages to model 
the light distribution and determine interesting fitted parameters. There seem to galaxy fitting 
parameters (at least) out there, GALFIT and GIM2D. I honestly cannot tell you which one is better. 

11.1 GIM2D 

GIM2D is a program that analyses galaxy images and spits out a galaxy-subtracted image (good to 
check the model with) and parameters. The use of SE is a crucial step in the process, 
website: 

\protect\vrule widthOpt\protect\href{http://www.hia-iha.nrc-cnrc.gc.ca/STAFF/Isd/g: 

11.2 GAEFIT 

GAEFIT also fits galaxy images. This was not initially intended to be used in automated batch 
analysis and SE is not a crucial step in the process, 
website: 

\protect\vrule widthOpt\protect\href{http://zwicky.as.arizona.edu\string'cyp/work/c 





Chapter 12 

SE use in the Literature 


In the literature SE is used for a myriad of problems. In this section the usage for detecting galaxies 
in the literature is discussed. It is not up to date completely, lots of people continue to use SE. Keep 
a look out in the literature. 


Casertano et al.l (l2000ar) used WEPC2 (U B V I) and a V+I detection image of the HDE-S with 


the following settings: THRESH 0.65 it,MIN_AREA 16 (0.05” pixelsize) and MIN_CONTAST 
0.03 . They used the weight maps made by drizzle for RMS background estimate. 


Brown et al.l (l2000ll used SE on photographic plates (UBRI) and made separate catalogs for 


each band. Classified as star if class > 0.7 in three bands or class > 0.75 in two bands or class > 
0.85 in one band. 

Rodighiero et al. ( 2000ll used SE on the HDE-N plus JHK band, detection was in the K-band 


excluding Vaucouleur-profiles. EEUX_RADIUS was used in photometry correction. 


Smail et al.l (1200 11) used SE on WEPC2 BVI with THRESH 1.5 cr 21.2 mag arcsec ^ 
MIN_AREA 10. Rejection criterion: k j 19 within 1-2 half light radii. 


and 


Kalirai et al.l(l200lh CEHT BVR data. Has a plot of magnitude -stellarity with a cut for the best 


star-galaxy separation. Mention of PSEex mentioned for correction of photometry. 

Williams ( 2001 ) run SE on Keckll data (BVI) to V = 26.5. All objects with Vj24 and CEASS 


i 0.2 is a galaxy. 

iRodighiero et al. ( 2001 ) run SE same as lEranceschini et'H^ ( 1998 ). SE is used to perform 
simultaneous slicing. Erom this, the PSE and the profiles for all objects are constructed. Galaxies 
are then selected as those objects resembling the Veaucouleurs profile convolved with the PSE... 
This is done for morphological classification. 

Trujillo et al. ( 200 lb used BVRI data from the NOT. I band detection with EOWTHRESHOED 


1.5(7 MINAREA 4. 2.5 kron radius photometry. 

Vaisanen et al. ( 2000b used SE on J,K data from the 1.2 m telescope. Used MAG_BEST and 


eyeball identification. 
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Zabludoff and Mulchaev ( 2000l) used SE in a standard config on rosat’s PSPC considering ev¬ 


ery thing with CLASS > 0.5 as a galaxy. Checked with isophotal area. MAG_BEST photometry. 

iKamhas et al. ( 12Q0Q ) used ground based data with MINAREA 5 (2”.3 pixelsize) and THRESH 
1 a. The SE catalog was then filtered for flagged objects. VLSB sample selection: /xq > 23 R a > 
3” isophotal area > 14 8 arcsec^. 


Ho gg et al.l (20001) used SE separately on UGR and K images, the EWHM of the smoothing 


Gauss was picked equal to the seeing. THRESH 1.2 cr and MIN_CONTAST 0.01 There was no 
Star/Galaxy separation. 


Volonteri et al.l (l2000ah used SE on the data of HDE (UBVI) south, in separate detections. 


The DEBLENDJVlINCON 0.01 DEBLENDJSTHRESH 32 smoothing Gauss 0.16 arcsec and 
MINAREA is equal to the seeing disk (!?!) THRESH E34. Photometry was ISOPHOT_CORR 
or aperture, depend i ng on d iameter ISO (1.2 arcsec limit). Stars haye class > 0.9 and I < 22 


Volonteri et al.l (I2000bl) used SE separately on the UVBI data of the HDE-S and later on a 


’meta-image’ of all filters combined, with THRESH 1.34 a and MIN_AREA 13 (seeing disk of 
arcsec^. Photometry was either ISOPHOT or aperture depending on isophotal diameter. The limit 
was 1.2 arcsec. To estimate the number of spurious sources, the image was inyerted and SE was 
run again on this image detecting only spurious sources. 

ISowards-Emmerd et al.l (l2000t) giye no details on their use of SE except that UB VRI detections 
we re done separate an d photometry was done on 20” apertures. 


Lubin et al. ( 2000l) used a image for the detection. This was made from three bands (BRi)using 


the Hale 5m tele s cope. 


Szalay et al.l (Il999bl) proposed the x^-images as the detection image. Seems only applicable 


wh en there is true multicolor i nform ation. 

iBroadhurst and Bouwens ( 2000 ) giye absolutely no details except that stars are badly fit by a 
red shifted galaxy spectra (Well DUH!) 

Castander and Lambl(ll999l)used S E to find candidates for yisual inspection. 


Gebhardt and Kissler-Patid (119991) used SE to search for Globular clusters in HST (VI) data. 


Globular clusters are those objects with an mag_err <0.1, V-I between 0 and 2, ellipticity < 0.5, 
EWHM between 1 and 4 pixels and ’nonstar’ classification. They estimate a constant contamination 
fro m background ga laxies. 


Poli et al.l (119991) use a ’meta-image’ from BVRI images. UBVRHK photometry was deter¬ 
mined from 1=23.25 isophotes and 2”.2 and 5” apertures. Half light radius, z and Magnitude 
rel ations. Not yery useful for galaxy id howeyer. 


Treu and Stiayellil (119991) used SE on NICMOS, WEPC and ground data. They distinguished 


between resolyed and unresolyed object by comparing with stellar EWHM’s. The limiting magni- 
tud e was the bin where all o bjects had a MAGERR_BEST < 0.15. 


Menanteau et al.l (119991) used SE on HST archiyal data for spheroidal g alaxies. They used bot h 


yisual identification and the concentration and asymmetry parameter from Abraham et al. llliS)- 
Th ey plot how to separate spirals f r om ellipticals using thes e parameters. 


Hashimoto and Qemlerl (119991) : iHashimoto et al.l (119981) used this concentration parameter to 


discern between early and late type galaxies. 
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Simard et al.l(ll999l) used Source Extractor on the Groth survey strip, THRESH E5 a, MIN_AREA 


10 and then fitting the li ght profiles of found galaxies. 


Nonino et al.l (119991) describe the SE program in detail, especially the detection process. The 


star galaxy separation is less stringent and said to favor complete star catalogs at class >0.5 
and more complete galaxy catalogs at class < 0.75. They present a flag list of SE. Parameters 
used:BACK_SIZE 64 (17”), THRESH 0.6 a. They use RMS images as weight maps (made by 
Weight watchers routine). MAG_AUTO is preferred and it is explained why. Good paper on the 
be havior of SE. _ 

van Dokkum et al. l2000i) touches only lightly on the use of SE but describes a CEEANing 


method to improve the resolution o f HSTAVEPC dat a. 


The study of high-z clusters 


ce et al. 


used EOCAS, SE and MDS ’find’ algorithm. 


SE is used on the K’ band images from the IRIM camera at the 4m Kitt peak. THRESH 1.5 u, 
{^K' = 22.2 arcsec^), 7x7 top-hat filter, MIN_AREA 35(0”. 15) The authors conclude that EOCAS 
an d SE do not show sig nificant differences in photometry, astrometry of classification). 


HiUcer et al.l (119991) present ground-based data and use SE in the V-band.. The background 


mesh is either 56x56 or 128x128 EWHM of the convolving gauss is either 1”.5 or 2” MIN_AREA 
5. The CEASS identification is taken to be accurate up to V <21. The limit was taken to be 0.35 
up to V = 22. Color is used to identify cluster membership. 


Smail et al.l (119981) naively used the MAG_BEST but used 1” diameter apertures on WEPC data 


and 3” on ground based data for galaxy colors. In the cases where Iband data was not available R 
ba nd data was converte d assuming R-I Ri 0.5 ± 0.2 


Teplitz et al.l (119981) used SE on NICMOS data with MIN_AREA 6 and used the EWHM of 


object, provided the were bright enough for a reliable estimate, to resolve stars from galaxies. All 
fai nter objects were assumed to be galaxies. 

Marleau and Simardl (119981) use SE for the initial detection of galaxies in the HDE and then de¬ 
termine a string of parameters with GIM2D to quantify the morphology. THRESH E5 (t, MIN_AREA 
30. DEBEEND_MINCON 0.001. 


Yan et al.l (119981) used NICMOS data EWHM 0.3” and a la detection and a \a analysis thresh¬ 


old. Plotofhalf-lightradius vs H magnitude with stars and galaxies. 

iBrown eraL ( 2001 ) analyze UBRI data from photographic plates from the SGP and UK Schmidt 
field. No details. 


Easano and Eilippil (119981) used SE on the HDE to find early type galaxies. THRESH 1.3 a. 
The selectionlimits for a galaxy were: Vqqq{STMAG) < 26.5 Npix > 200 CLASS < 0.6. 
lEasano et al.l ( 1998 ) used this data to establish a SB- effective radius relation for the early type 
galaxies. Note that they chose to use the V band for identification due to better S/N. 

? used their own program to make a survey of southern hemisphere faint galaxies. More 
information than surmisable here. 


Hogg et al.l (I1997D use SE on U and R band images. No specifics are given except that the 


det ection was done on images sm oothed with the PSE. 


Bertin and Dennefeldl (I1997D used SE on digitized Schmidt plates. It contains a good descrip¬ 


tion of the usage of SE on photographic plates and the particular problems encountered. 
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Gardner et alJ (Il996r) use SE on BVIand K band ground based images. They use SE on all of 


them and note that I-K and B-I color are good separators for star galaxy for all except the bluest 
obj ects. All identified ga laxies up to certain limits were confirmed by eye. 

Lanzetta et al. ( 1996ll use SE on HST/HDE images. Detection is in the I band. The EWHM 


smoothing was 0.12 arcsec (approx width PSE) to aid detection of faint sources. Einally there 
is a re ference as to why you should bother to smooth again with approximately the PSE: llrwin 
( 1985li .MIN_AREA 10 THRESH 1.4 a and surprisingly a CLEAN_PARAM pf 2.0 to get rid of 
closely-packed objects (taken to be one object) BACK_SIZE was an unusual 41x41 pixels. As 
detection was done in the I band image, photometry was done in dual mode for the rest of the 
bands. 


Small et al.l (Il997h produced catalogs on 10 clusters images by HST. Detections as in the red¬ 


dest band available (either E702W o r E814W). The modifie d SE code to produce the concentration 
index introduced hv lAhraham et al.Nl994h : lAhraham et akl and the contrast index which measures 
the fraction of light in the brightest 30% of its pixels. Object were however checked visually. 
THRESH 1.3 cj is the only parameter mentioned. 





























Chapter 13 

Known Bugs and Features 


Figure 13.1: Comic taken from http://www.phdcomics.com/ No royalties paid. Hope that’s ok. 
Please don’t hurt me. 

SE was, still is, always will be in develoment. After all it’s pretty much one guy coding it all up 
(Emmanuel Bartin if you’ve missed that). So there are still outstanding issues. Here are the ones I 
found. They’ll be mentioned throughout but here’s the list: 

1. The PSE fitting routine does not work (yet?). 

2. Neither does the galfitting routine. But hey there are follow-up scripts for those. 

3. The EEEIPTICITY and all the moments is determined from the smoothed image. This is a 
faeture not a bug but it is not what people intuitively expect SE to do. 

4. The first pixel’s coordinates is 1,1. If you’re thinking in terms of matrices, this can be an¬ 
noying. 

5. The documentation will always lag the program (but to claim this is unique to SE...) 

6. Output parameters asked for are important. This is an SE feature. If no ASSOC output 
parameter is asked for, the assoc module is not run. If no Kron-dependent parameter is in the 
param file, this aperture is not drawn on the APERTURE image. It’s a feature you need to 
be aware of. 





Chapter 14 

SE parameter additions 


It’s also discussed in the manual but there is a possibility that the user can define his/her own output 
parameters in SE. In order to do that, you need to do the following things: First ask yourself if you 
can’t do it with ANY of the parameters above or combinations thereof. After all there is a lot of 
them and it a lesser headache. 

If not, then define what this parameter should be and then modify the analyze.c, types.h and param.h 
as follows: Define your parameter after the definition of the FWHM function in analyse.c (after 
line 371 in version 2.2.2). Define the parameter name in types.h and param.h like so: 

In types.h 

float cone; /* IRS concentration index */ 


and in param.h: 

{"CONCENTRATION", "Abraham concentration parameter", 

&outobj.cone, H_FLOAT, T_FLOAT, "%8.3f",""}. 

This way, if you include CONCENTRATION in the file given to PARAMETERS _NAME, then 
the concentration will be calculated and put in the value ’cone’. 


This parameter is the concentration param eter as defined inlAbraham et ah ( 1994 ); Abraham et al 


and previously implemented by Ian Small in dCouch et al. ( 1998ll . I’ve modified it so that it works 
in version 2.2.2. 

The contrast and asymmetry parameters are given as examples later on. 


/* Abrahams concentration index calculated - IRS/BWH */ 

/* fraction of light in central 30% of the objects area, measured 
in an ellipse aligned with the object and having the same axis ratio */ 
if (FLAG(obj.cone)) 

{ 
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double xm,ym,dx,dy,rv,cv,amp; # local variables 
double AA,BB,CC,DD,EE,FF,sintheta,costheta ; 

xm = obj->mx; # get the median x position from the object struct 
ym = obj->my; # get the median y position from the object struct 
amp = tv/{2*PI*obj->a*obj->b*obj->abcor); 
rv = cv = 0.0; 

for {pixt=pixel+obj->firstpix;pixt>=pixel;pixt=pixel+PLIST(pixt,nextpix)) 
# run through all the pixels in the object struct 
{ 

dx = PLIST(pixt,x) - xm; # calculate dx and dy 
dy = PLIST(pixt,y) - ym; 

pix = PLIST(pixt,value)<prefs.satur_level? 

PLIST(pixt,value) 

: amp * exp(-0.5*(obj->cxx*dx*dx + obj->cyy*dy*dy + obj->cxy*dx*dy) 
/obj->abcor); 

costheta= cos(obj->theta*PI/180 . ) ; 
sintheta= sin(obj->theta*PI/180. ) ; 

EE = (obj->b>0.0)?sqrt(obj->npix*obj->a/{PI*obj->b) ) :1.0; 

FF = (obj->a>0.0)?sqrt(obj->npix*obj->b/{PI*obj->a)):1.0; 

AA= costheta*costheta/EE + sintheta*sintheta/FF ; 

BB= CC= sintheta*costheta* (1.0/EE - 1.0/FF) ; 

DD= + sintheta*sintheta/EE + costheta*costheta/FF ; 

if {{AA*dx+BB*dy)*{AA*dx+BB*dy)+{CC*dx+DD*dy)*{CC*dx+DD*dy) < 0.09) 
cv += pix; 

if {{AA*dx+BB*dy)*{AA*dx+BB*dy)+{CC*dx+DD*dy)*{CC*dx+DD*dy) < 1.00) 
rv += pix; 

} 

obj->conc = (rv>0.0)? (cv/rv): 99.0; 

} 


Chapter 15 

Examples 


In this chapter I intend to show the workings of Source Extractor using the Hubble Deep Field and 
Hubble Ultra Deep Field. Of course this does not cover the many many uses of SE but it should 
give you an idea. 

15.1 The Data: the Hubble Deep fields 

The Hubble deep fields were asfronomy’s firsl look af fhe high-redshiff universe. As a resulf fhese 
fields have been analysed fo deafh and are well known. These also saw some of fhe firsl quick-look 
papers lhal almosl always used SE for fhe slalislics. Also, high-end science producls are freely 
available. 

The Hubble Deep Field (Norlh) can be found here: 


\protect\vrule widthOpt\protect\href{http://www.stsci.edu/ftp/science/hdf/hdf.html' 
\protect\vrule widthOpt\protect\href{http://stdatu.stsci.edu/hst/hdf/v2/mosaics/x4( 

I am using fhe files which can be found here: 


The Hubble Ullra Deep field is a Advanced Camera for Surveys producl. This is a much bigger 
field and hence file. I’ll use fhese lo illuslrale fhe capabilily of SE lo habdle larger files^ 

\protect\vrule widthOpt\protect\href{http://www.stsci.edu/hst/udf}{http://www.stsc: 
\protect\vrule widthOpt\protect\href{http://archive.stsci.edu/prepds/udf/udf_hlsp.1 

*I realise that computers are getting bigger and better so the definition of “a large file” is constantly changing but the 
UDF is freely available and easy to obtain. 
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15.2 Running the SE default 

Simpelest place to start. Put the default files in the config directory from source extractor in the 
same directory as the f814_mosaic_blk4.fits file. Run SE wifh: 

> sex f814_mosaic_blk4.fits 

And fhaf is if! This has fafhered^ a cafalog named ‘fesf.caf’. 

Since fhe defaulf.param oufpuf was generafed fhis confains fhe following parameters: 


# 

1 NUMBER 

Running object number 


# 

2 ELUXERR_ISO 

RMS error for isophotal flux 

[ count] 

# 

3 ELUX_AUTO 

Elux within a Kron-like 


elliptical aperture 

[count] 


# 

4 ELUXERR_AUTO 

RMS error for AUTO flux 

[ count] 

# 

5 X_IMAGE 

Object position along x 

[pixel 

# 

6 Y_IMAGE 

Object position along y 

[pixel 

# 

7 ELAGS 

Extraction flags 



And 1714 enfries. The checkimages can be obfained as well by running fhis: 

sex f814_mosaic_blk4.fits -checkimage_type APERTURES,BACKGROUND 
-checkimage_name aper.fits,back.fits 

fhese are in figure (TSTI 

Figure 15.1: The Hubble Deep Eield North (E814W) and the various checkimages when SE is run 
with the default settings. The apertures file shows the KRON apertures around the objects. 


15.3 Stars and Galaxies 

Now we go into the default.param file and comment out everything except the MAG_AUTO and 
the CEASS_STAR parameters. If we plot these two values (still using default values for the source 
extraction) we get figure ??. 

Figure 15.2: The MAG_AUTO and CEASS_STAR parameters of the Hubble Deep Field. The 
MAG_APER has strange values as the zeropoint setting is wrong. Most of the objects are classified 
as galaxies but the classification becomes random for the dimmer objects. 

Note how similar it looks to figure 18.31 What is you need to take away is that the star/galaxy 
separator is okay to give an indication but that it becomes a random number genereator pretty quick. 

^sex, fathered...get it? Oh the puns...ugh. 
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There are not that many stars in the HDF yet there are quite a number of objeets elassified as sueh. 
As with most elassifieation sehemes, SE finds out if an objeet is extended or not. Not extended, 
must be a star... 

15.4 Photometry 

15.5 Typieal Radii, souree sizes 
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Appendix A 

’Drizzle’ and RMS weight images 


Suppose you’re a HSTAVFPc user and you want to use SE on the data. The program ’drizzle’ which 
combines a series of exposures produces also weight maps of the ’drizzling’. These weight maps 
can be used as MAP_RMS maps if processed as follows: weightmap = sqrt[drizzieteight^ap) wi* 
Fa a certain correction parameter to get pure, uncorrelated RMS noise. This value depends on the 
area over which the noise is determined. Casertano et alJ ( 200Qh ) estimated this in their appendix 
and came to this: If the pixfrac of drizzle is p and the scale of the output pixels is s (0.5 half of the 
original pixels) then the ratio between uncorrelated noise and the drizzle map value is: 

S 1 5 

sqrt{FA) = (-(1 - --) , if s <p 

p 6p 

sqrt{FA) = 1 - ^s/p) , if s> p 


If you want to use this trick, please go over this bit in ICasertano et al. (2000b) to be sure. 

The resulting map can be used with WEIGH_TYPE MAP_RMS and thus takes in account all 
the weird thing drizzle did to the data (per pixel different contributions from pixels with maybe 
different exposures). 














Appendix B 

SE parameters 


Below the settings of SE are listed with their function and the best way to use them in a crowded 
field. These values were obtained after using SE on a series of simulations of Sextans A (north 
pointing wf2) and the HDE. The success rate was measured by the number of objects classified as 
galaxies in fhe simulations acfually found in fhe HDE wifh fhe same sellings. 


CATALOG_TYPE ASCII # This means no header with the values 

# in different colums is printed. 

PARAMETERS_NAME /net/bartoli/bartoli/data2/Programs/default2.param 

# file with the names of the parameters to be 

# printed in the catalogs. 

DETECT_TYPE CCD # type of image (alt PHOTO) 

DETECT_MINAREA 16 # Minimum area of connected pixels in an object. 

EILTER Y 

FILTER_NAME /net/bartoli/bartoli/data2/Programs/sextractor2.1.6/config/gauss_4.0_7 

DEBLEND_NTHRESH 32 # the number of thresholds the intensity 

# range is devided up in. 32 is the most 

# common number. 

DEBLEND_MINCONT 0.005 # percentage of flux a separate peak must 

# have to be considered a separate object. 

CLEAN Y # Should objects close to bright object be 

# removed from the catalogs? 

CLEAN_PARAM 1.0 # measure for 'clean'ing. wish I knew what it ment... 
PHOT_APERTURES 5,10,20,30 # Fixed apertures in pixels 
PHOT_AUTOPARAMS 2.5, 3.5 

SATUR_LEVEL 200000. # saturation occurs at? 
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MAG_GAMMA 4.0 # Emulsion response slope 

# SE doesn't work without this for some reason 
GAIN 7.0 # number of photons / counts ratio 

# if run over a counts per second image then the 

# gain*exp time is used. 

PIXEL_SCALE 0.05 # arcsec 
SEEING_EWHM .17 # 

STARNNW_NAME /net/bartoli/bartoli/data2/Programs/sextractor2.1.6/config/defai 
BACK_SIZE 32 # size in pixels of the area used to estimate the 

# background 

BACK_FILTERSIZE 1 # type of filter used in the background 
BACKPHOTO_TYPE LOCAL # local estimate or an estimate for the entire 

# image at once... 

BACKPHOTO_THICK 32 # 

#CHECKIMAGE_TYPE APERTURES # output image of SE. 

#CHECKIMAGE_NAME Roza_F555W_wf3.fits # name of that image. 

CHECKIMAGE_TYPE NONE 

MEMORY_OBJSTACK 10000 # memory parameters 
MEMORY_PIXSTACK 1500000 # defaults work best 
MEMORY_BUFSIZE 1024 
#SCAN_ISOAPRATIO 0.6 
VERBOSE TYPE NORMAL 
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